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SUMMARY 



The application of plasma chemistry for removal of Volatile Organic Compounds (VOC) 
in the forest products industry is studied in this work. An extensive research was done to 
understand the physical and chemical processes with VOC in pulsed corona and gliding arc 
discharges at conditions typical for the exhaust streams in paper mills and strandboard press 
vents. 

A linear stability analysis of the low-current gliding arc discharge in the Uransitional 
regime was performed to explain the observed experimental effects in VOC removal. It was 
shown that the gliding arc remains stable during the evolution and gradually transforms into a 
more non-equilibrium one. The low-ciurrent arc discharge can propagate with the effect of 
"overshooting" at which the gliding arc extinguishes long after its maximum power has been 
achieved. Analytical and numerical solutions explain the general behavior of the low-current 
gliding arc and are in a good agreement with our experiment. 

A semi-empirical approach of coupling physical parameters of a pulsed corona discharge 
with its chemistry is proposed. A fraction of high-energy electrons produced by the discharge was 
found using reaction mechanism of NOx formation in pure dry air. An assembled plasma- 
chemical elementary-reaction mechanism for methanol and dimethyl sulfide showed the validity 
of this approach. A Projected Energy Cost (PEC) was proposed to generalize the obtained 
experimental data. The optimum PEC for VOC removal exists at 2% of absolute humidity. The 
worst removal rate was observed at 130°C. The results were used for the process improvement in 
the mobile pulsed-corona pilot plant that is currently under construction. A concept of the plant is 
discussed. 

X 



1. INTRODUCTION 



1.1. Current Technoloaes for VOC Abatement 

The forest product industry is coming under increasing pressure from the U.S. 
Environmental Protection Agency as well as state and local agencies to control emissions from 
pulp mills and wood products plants. Incineration (and other high temperature processes) is 
widely used environmental-control technologies for concentrated volatile organic compound 
(VOC) emissions; however they may not be appropriate for the dilute VOC process emissions 
commonly found in the forest products industry. Although direct incineration of dilute VOC 
streams is technically feasible, it will typically be much more expensive than low-temperature 
technologies that avoid heating the entire process stream with supplemental fiiels, especially in 
the case of high-volume streams. Potential low-temperature purification technologies include 
low-temperature catalysis, biofiltration or biotreatment, and low-temperature, non-equilibrium 
plasmas. Low-temperature catalytic purification processes have disadvantage of catalyst that are 
subject to plugging and poisoning. Thus, these systems tend to have higher maintenance and 
replacement costs than the plasma technologies. Biofiltration requires high surface areas and long 
residence times for VOC capture and destmction. This results in substantially larger process 
volimies than those that can be handled by plasma technologies. Biofiltration may also require 
more carefiil process control in order to maintain bioactivity for VOC control. Biological 
treatment is a substance-sensitive. In this thesis the biofiltration will not be considered for 
comparison. 



I 



2 

Low-temperature, non-equilibrium plasmas (also called "non-thermal" plasmas) are an 
emerging technology for abating diluting VOC emissions. These plasmas may be produced by a 
variety of electrical discharges or electron beams. The comparison of non-thermal plasma 
technologies with incineration techniques is given in Table I. 

The basic feature of plasma technologies is that they produce plasma in which the 
majority of the electric energy (more than 99%) goes into production of energetic electrons, 
instead of heating the entire gas stream. These energetic electrons produce excited species, free 
radicals and ions) as well as additional electrons through the electron impact ionization, 
excitation and ionization of the background molecules. These excited species, in tura, oxidize, 
reduce, or decompose the pollutant molecules. This is in contrast to the mechanism involved in 
thermal processes (such as plasma torches or fiimaces. Regenerative Thermal Oxidation (RTO) 
and several chemical techniques) which require heating the entire gas stream in order to destroy 
pollutants. In additional, the low-temperature plasma technology is highly selective and has 
relatively low maintenance requirements. Its high selectivity results in relatively low energy costs 
for emissions control while low maintenance keeps annual operating expenses low. Furthermore, 
these plasma discharges are very uniform and homogeneous (except gliding arcs) which results in 
high process productivity. 

Specific niches of incineration and plasma technologies are clearly seen in Table 1. As 
mentioned above, RTO is more effective at high concentrations of VOC. Plasma is normally 
used al lower concentrations of VOC and in the cases where combustion is not effective, such as 
removal ofNOx or sulfiir-containing compounds. Scrubbing of exhaust gases with water might be 
required for both techniques, but some plasma reactors can be designed as a single unit with a 
scrubber to reduce the equipment cost. Moreover, it was shown in this thesis, that higher than 



1000 ppm concentrations of VOC can be successfully treated in pulsed corona reactors. In the 
foreseen future, due to new environmental regulations plasma may take a substantial share of the 
gas abatement market. 



TABLE I. COMPARISON OF RTO AND NON-THERMAL PLASMA TECHNOLOGIES 



\3us ircsnneiic 
method 


Regenerative Thermal 
Oxidation (RTO) 


Non-equilibrium Plasma Technologies 


Applications 


Gas Streams with high 
concentration of VOC 
(> 1000 ppm) 


Gas Streams with low concentration of 
VOC (< 1000 ppm); removal of NOx 
andSOj 


\Iethod Selectivity 
or Electrical 
Efficiency 


Not selective, heating of 
the entire flow up to high 
temperatures (> 1000 "C) 


Up to 99% of electrical energy could go 
to production of active species, almost 
no heating of gas 


Safety 


Potential problems of major 
accidents in case of lack of 
treated air flow, or control 
power outage. High 
temperatures. 


Safe, but employs high voltage 
equipment 


Equipment and 
Maintenance Cost 


High 


Mediimi 


Need for Additional 
Gas/Liquid Supply 


Needs considerable 
amounts of natural gas to 
support combustion 


Might need some scrubbing water 


Potential Problems 
with Byproducts 


Sulfur-containing 
compounds must be treated 
additionally 


Dangerous by-products may be present 
in the stream; process is sensitive to the 
kind of VOC; complex compounds 
require longer residence time 
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1.2. Overview of Plasma Technologies for Abatement of VOC 

Low-temperature, non-equilibrium plasma processes have been shown to be effective in 
treating a wide range of emissions including aliphatic hydrocarbons, chlorofluorocarbons, methyl 
cyanide, phosgene, formaldehyde, as well as sulfur and organophosphorus compounds. Although 
the products of these plasma processes are virtually undistinguishable from incineration products 
(COj, H2O, SO2, etc.), chemical reactions occurring in these technologies are substantially 
different than in incineration ([1 1], [62]). 

Large electric fields in discharge reactors create conditions for electric breakdowns during 
which many electron-ion pairs are formed. Primary electrons are accelerated by the electric field 
and produce secondary ionization, and so forth. Created excited species, atoms, radicals, 
molecular and atomic ions, electrons and radicals are capable of interacting with VOCs to a 
certain degree. Some processes are more effective than the others. The most important are those 
leading to formation of OH radicals. Typically, the major mechanisms controlling the chemistry 
of VOCs fall in five major categories. 

1.3. Radical-Indnced VOC Dcstrucrion 

Hydroxyl radicals are formed through a three-step mechanism initiated by energetic 
electrons, which ionize any specie (M): 

e + M^VT + e + e 
and followed by the fast charge transfer 

M* + H20^M + H20* 
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The positive water ions react rapidly with neutral water molecules to form the following 
radicals: 

H20^ + H2O + OK 

This mechanism is especially important in high voltage discharges, such as pulsed corona 
and dielectric-barrier discharges, and electron-beam systems. 

Excited and ground state atomic oxygen (0('D) and 0(^P)) species are also important 
radicals for producing hydroxyl radicals 

0('D) + H20->'0H + '0H 

These atomic oxygen species also can produce ozone O3, which is particularly active in 
liquid-phase oxidation. 

1.4. Direct Electron-Induced VOC Destruction 

This mechanism usually only occurs in strongly electronegative gases, when the electron 
affinity of a gas specie is comparable to its dissociation energy. For example: 
e + CHCl3->'CHCl2+Cr 

or 

e + CCl4->'CCl3 + Cr 

Usually a negatively charged particle is lost in fast ion-ion recombinations. Higher 
electric fields favor such type of mechanism. 
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1.5. Direct Ion-Induced VOC Destruction 

The direct ion-induced decomposition of VOCs is similar to the hydroxyl-formation 
mechanism discussed earlier. When either the VOC species or their intermediate products 
produced by another destruction mechanism have a low ionization potential, the charge exchange 
process can promote further decomposition reactions. 

Direct electron impact ionization of nitrogen molecules occurs in collision with an 
additional nitrogen molecule: 

N2 + e Ni* + e + e 

If the ionization potential of a hydrocarbon radical (produced in an earlier VOC- 
destruction reaction) is very low, the charge exchange process can provide dissociation of organic 
molecules as well. Methanol decomposition by the positive nitrogen ion is a good example of this 
direct ion-induced mechanism by the reactive charge exchange process: 

N2* + CH3OH -> CH3* + 'OH + N2 

For all three abovementioned mechanisms of VOC destruction, energy consumption for 
the gas purification process is directly related to ionization energy. The higher the electric field, 
the higher the average electron energy, and the greater than the fraction of the total discharge 
power that will go into ionization, increasing the efficiency of the VOC-destruction process. This 
fact is one of the key points for applying cold discharge technology to environmental control 
processes. 
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1.6. Water Droplet- and Cluster-Enhanced VOC Destruction 

The effect of humidity on plasma-chemical efficiency was studied by several authors. 
While several reports have noted that the humidity has little or no effect, the other reports say that 
the effect of humidity is detrimental for the removal efficiency. In this thesis the independent 
investigation of the role of water was carried out. 

Ions produced in an atmospheric-pressure discharge are natural nuclei and would 
stimulate water condensation or cluster formation. The effect of droplets is also dual. Liquid 
phase catches positive ions and ozone molecules. Captured species induce formation of OH 
radicals in droplets. The radicals very efficiently destroy VOC molecules when they encounter 
the droplet. 

The experimental data presented in this thesis showed some interesting patterns of 
dependence of VOC destruction and removal efficiency on water concentration. It has to be 
noted, however, that even though the diverse role of water for VOC treatment is recognized, 
quantitative models of all these effects are yet to be developed. 

1.7. Ultraviolet VOC Destruction 

Ultraviolet (UV) radiation in plasma can also be effective in selective dissociation of 
VOC (Rusanov and Fridman 1984). Although this effect is not a primary decomposition 
mechanism in non-equilibrium plasmas, plays a very important role as an ancillary VOC- 
destruction process. Typically, LTV radiation emitted by the plasma discharge breaks molecular 
bonds. 
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While multiple VOC-destruction mechanisms take place in plasma processes, the primary 
mechanism involves the dehydrogenation of an organic site by a hydroxy! radical (OH) followed 
by attachment of oxygen as described in the following sequence of reactions. 

R-H-HsCH + OH" R-H-H2C+H2O 

R-H.H2C+O2 R-H-H2C-O-O" 

The organoperoxide radical formed in the second reaction then reacts with another 
organic specie to form an organoperoxide while dehydrogenating the second organic molecule. 
R.H-H2C - O-O* - R-H-H2CH R-H-H2C - O-OH + R-H-H2C • . 

The organoperoxide is unstable and decomposes to a new specie of one less carbon atom 
than the original compound and either formaldehyde or formaldehyde radical, which is rapidly 
oxidized by oxygen to carbon dioxide and water. 

R-H-H2C - O-OH -> R'" '^-H + (H2CO or HCO" ) 

H2COorHCO" +0i^ CO2+H1O 

Aromatic and unsaturated organics follow a slightly different decomposition mechanism 
but are both eventually decompose to carbon dioxide and water. In general, the unsaturated 
compounds are easier to decompose. The chain of plasma-chemical reactions shown above can 
be also initiated by atomic hydrogen or nitrogen, but these branches are significantly less 
probable. 

It is seen that presence of considerable amount of water in the treated air is essential for 
effective destruction of VOC. The electron energy is the primary factor in efficiency of all major 
plasma-chemical mechanisms of VOC destruction. From this point of view the quality of plasma 



produced by electrical discharges or electron beams is different. Table II summarizes the features 
of major plasma technologies for VOC abatement applications. 



TABLE II. COMPARISON OF NON-THERMAL PLASMA TECHNOLOGIES ([!]) 



Plasma 
Technology 


Average 
Electron 
Energy 


Capital 
Cost, 
SAV 


Unit 
Capacity, 
kW 


Energy Con- 
sumption, 
kWh/m' 


Known 
Applications 


Technical 
Limitations 


Electron 
Beam 

Technology 


5-6 eV 


2.0 


1,000 


10-70 


Removal of 
SOi, NOx in 
thermal power 
plants 


Initial energy of 
electrons in the beam 
is limited by 1.5 MeV 
due to x-ray 
production 


Gas-Phase 
Corona 
Reactor 
(GPCR) 


2-3 eV 


2.1 


20 


6-70 


Electrostatic 
Precipitators 


The voltage cannot be 
as high as in pulsed 
corona due to spark 
formation; packed 
bed limits flow rates. 


Pulsed 
Corona 
Discharge 


3-5 eV 


1.0 


10 


20-150 


Removal of 
polyatomic 
hydrocarbons 
in aluminum 
plant 


Voltage is limited by 
the electrical network 
coupling and 
capacitance and 
inductance of the 
plasma reactor 


Dielectric- 
Barrier 
Discharge 
(DBD) 


2-3 eV 


0.2 


0.1 


20-200 


Removal of 

S02,N0x, 

industrial 

ozone 

generators 


The voltage cannot be 
as high as in pulsed 
corona due to spark 
formation; low unit 
capacity 


Gliding Arc 


1-1.5 eV 


O.l 


120 


20 - 250 


Laboratory 
study on 
removal of 
methanol, o- 
xylene 


Discharge is non- 
uniform, high 
temperatures of the 
gas, possible NOx 
generation. 
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Electron beams produce plasma of the best quality with average electron energy of 5-6 
eV, Pulsed corona discharges stand very close (3-S eV). But corona discharges have advantage of 
lower capital cost. Generally, both technologies have aheady found applications for gas treatment 
form NOx and SO2 emissions in coal power plants. Capability of treating large gas streams make 
e-beams and pulsed corona discharges the primary candidates considered for VOC abatement. 

1.8. Goals of the Forest Project 




The goal of this project is to test four candidate VOC control technologies based on non- 
thermal plasma, and demonstrate the best technology on a mill. They are: 1) dielectric packed bed 
corona discharge (GCPR); 2) electron beam discharge; 3) non-thermal gliding arc discharge; and 
4) heterogeneous pulsed corona discharge. The development of the last two technologies (gliding 
arc and pulsed corona discharge) will be discussed in this thesis. 

The target streams for this study are the vent streams from brownstock washers and 
oriented strandboard presses and dryers. Very little data exists on contaminants of interests in the 
forest products industry and test conditions were far from typical industrial streams. 

Press vents are especially attractive for application of non-thermal plasma due to it 
intermittent nature. Most of the time press vents operate with very little concentration VOC. 
However, a spike of VOC is emitted when the press is open. This spike can create problems for 
thermal control technologies. They are wasting the energy during a great fraction of the operation 
cycle. Non-thermal plasmas, on the other hand, are essentially "instant on". This opens the 
opportunity for operating non-thermal plasma device "off" or "idle" until just before the spike 
arrives and then turn on (or increase) the power for the duration of the spike. Significant energy 
savings are expected from this type of "pulsed" operation. 
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The laboratory experiments will determine which of the four plasma technologies would 
be most suitable for treating specific VOC emissions of concern to the wood products industry. 
Simulated VOC emissions (typical of current operations) will be prepared in the laboratory by 
blending a combination of gases and liquids to obtain realistic gas compositions, humidity, and 
temperatures as shown in Table III. Concentrations during spikes are described as the worst-case 
scenario. The feed gas compositions will be confirmed chromatographically. 



TABLE III. TYPICAL STREAM COMPOSITIONS CONSIDERED IN THE PROJECT 





HVLC 
Brownstock 
Washer Went 
Emissions 


Oriented 
Strandboard Press 
Vent Emissions 


Oriented 
Strandboard Dryer 
Vent Emissions 


Gas Composition: 








Acetaldehyde 






7 ppm 


Dimethyl Sulfide 


2 ppm 






Methanol 


83 ppm 


25 ppm 


8 ppm 


Acetone 


3 ppm 


1 ppm 


28 ppm 


Terpenes 


209 ppm 


4 ppm 


16 ppm 


1,2,4- 

Trichlorobenzene 






2 ppm 


Acrolein 






2 ppm 


Humidity 


100%RH@110°F 


70%RH@ lOOT 


90%RH@200'*F 


Process Conditions: 








Temperature 


110°F 


lOO'l^ 


200°F 
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The demonstration of the selected technologies will be conducted at mills selected by 
Georgia-Pacific Corporation. These demonstrations will allow us to validate performance over 
typical mill operating conditions, determine equipment performance over an extended (2-4 
months) period, and determine maintenance requirements, which will govern the true cost of 
ownership. 

1.9. Development of the Trailer-Mounted Pilot Plant 

The ultimate goal of this project is testing the prospective pulsed corona technology at 
actual on-site operating conditions of the process plants in paper industry. A Louisiana paper mill 
plant of Georgia-Pacific Corporation in Port Hudson was selected for the testing, scheduled for 
Fall 2002. The testing on streams from the Brownstock Washer Vents will be the introductory 
phase of the pilot testing. In case of successful operation in Port Hudson, the technology should 
be demonstrated on exhaust streams from the Oriented Strandboard Press Vents and from the 
Oriented Strandboard Dryer Vents. The major compounds being treated in all three cases are 
essentially the same. However, their concentrations, temperatures and absolute humidity of the 
exhaust gases might differ. 

Recognizing the need to evaluate nonthennal plasma technology for application in the 
forest products industry. National Council for Air and Stream Improvement (NCASI) reviewed 
the technology options. 

Table IV shows the input data on VOC used in the report [1] issued by NCASI in 1999. 
The expected gas composition of in the paper mill considered in the report appeared to be 
significantly below the real concentration of noxious species in the vent stream ("the worst-case 
scenario"). 
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TABLE IV. COMPARISON OF DATA ON BROWNSTOCK WASHER VENT EMISSIONS 



Composition 


NCASI 


"The worst case scenario" 


Dimethyl 
Disulfide 


2 ppm 


20 ppm 


Dimethyl sulfide 




1 727 ppm 


Methanol 


83 ppm 


2330 ppm 


Acetone 


3 ppm 




Terpenes 


209 ppm 


62 ppm 


H2S 




Up to 100 ppm 


Conditions 


Temperature 


103 F 


150 F 


Relative Humidity 


100% 


100 % 



The following differences of "the worst case scenario" must be underlined because they 
constitute additional challenges for the successful pilot plant operation. 

Concentrations of methanol and dimethyl sulfide exceed the expected value about 20 
times. Such concentrations of VOC are very challenging for the corona unit to be designed. 
Typically, corona is most feasible for concentration of VOC in the range of 100 ppm. 

Methanol, presented in such amounts, must be removed in the earliest stage. Otherwise, it 
will degrade the performance of the corona by consuming a major fraction of the active radicals 
produced by corona. 
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The stream contains dangerous amounts of HjS. It is about 100 times larger than the 
threshold limit imposed by regulations. Additional technical measures would be required to 
maintain the safety of the pilot plant operation. 

Presence of large amount of water in the stream makes it difficult to predict performance 
of the pilot plant at high absolute humidity and at elevated temperatures. Water molecules change 
the property of the discharge and also might take effect on chemistry. 

The evident lack of experimental data for the given substances and on performance of the 
pulsed corona at high absolute humidity, high concentration of VOC, and at elevated 
temperatures, motivated the research presented in this thesis. 

Actual demand of capacity of the process gas treatment equipment in the paper mills 
exceeds 100,000 standard cubic feet per minute. It is expected that such treatment facility will be 
built from basic units of much smaller capacity (about 1:100). From the purpose of technology 
demonstration, it is sufHcient to test in field just this unit for the time comparable with actual 
lifetime of the equipment or with maintenance cycle. It was proposed to test the pilot plant for at 
least 3 weeks in the regime of continuous operation. Table V sunmiarizes the chosen design 
parameters of the pulsed corona pilot plant. Research and development of the pilot plant was 
done in collaboration of several organizations. 

University of Illinois at Chicago (later Drexel University) leaded the effort. The author 
of this thesis supervised the design of the pilot plant hardware and did the engineering work on 
integration of the pilot into the trailer, proposed the program of on-site testing. Research on 
physics and chemistry of pulsed corona discharge was done in collaboration with Mario Sobacchi 
and Dr. Andrey Romanets, and under the supervision of Professor Alexei Saveliev. The 
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feasibility study and theory development of the gliding arc discharge was completed in 
collaboration with Mario Sobacchi and Inga Kuznetsova. 



TABLE V. OPERATIONAL PARAMETERS OF THE PULSED CORONA PILOT PLANT 



Gas flow rate 


Up to 600 m^/hour 


Consumed Electric Power 


Up to 10 kW 


Consumption of water for the scrubber and 
corona units 


Up to 0.7 m-'/hour 


Consumption of compressed air for the spray 


40 m^/hour 


Pressure of air inside the pilot plant 


6 negative inches of water (vacuum). 



EGOS Ltd (Kurchatov Institute of Atomic Energy, Russia) developed a prototype corona 
ceil that was used to build a pilot plant. They built the pilot plant hardware. 

Georgia-Pacific Corporation paper mill plant in Port Hudson was chosen as a place for 
the testing. 

Argonne National Laboratory provided the expertise on project management and 
interfacing with local engineering teams. 

Pacific Nortiiwest National Laboratory developed a packed-bed corona pilot 
installation for testing according to phases II and III of the project. 

Joint efforts of the team succeeded in conceptual 3D design of the mobile pilot plant 
(Figure 1) and in development of logic of its operation. 
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The pilot plant will be mounted on a 48-foot long single-drop flatbed trailer to minimize 
the overall height. A standard 10-foot long mobile office will serve as a control room and an 
analytical laboratory. 

As seen in Figure 1, the pilot plant seems very tall and narrow. Its center of gravity is 
high. Unfortunately, the vertical dimension is the least customizable. It is set for the given mass 
transfer rate in the scrubber and in corona, and is less than the height of bridges on the highways. 
To avoid flipping of the trailer in the sharp turns, it must be much heavier than the plant. For the 
selected trailer this condition holds. 

The width of the installation should not exceed the width of the trailer and should allow 
some space around the plant for service work. Narrow cargo helps stabilize the trailer at high 
speeds on highways and at side winds. 

The installation is centered with respect to the longitudinal axis of the trailer. For this, the 
preliminary design of the pilot plant was modified to move very heavy transformers from the side 
to the fh)nt of the Urailer. 

1.9.1. Pescriptioii of the Design of the Mobile Pul sed Corona Plasma Pilot Plant 

According to the concept of operation, the pilot plant incorporates three different cleaning 
technologies, represented by individual modules: a) scrubber; b) pulsed corona modide; c) 
activated carbon bed/mist separator module. All three modules are assembled as a single unit 
with a common enclosure. Auxiliary equipment, including a waste water pump, air blower and 
high voltage transformers, will be mounted outside the enclosure directly on the flatbed. 
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Scrubber 

The air with VOC is supplied for treatment to an incineration facility by a major vent 
pipeline. To avoid potential hazards due to leakage, the air in the pipe is always under slight 
negative pressure (vacuum) of 6 in of water. The pilot plant works on the same principle. The 
flow through the modules is sustained by a powerful auxiliary blower capable of producing up to 
more than 10 negative inches of pressure with the required flow rate. It is located downstream of 
the pilot plant on the exhaust pipe. A 6-in diameter flexible hose will connect the pilot and the 
vent. If the blower should fail or the hose should leak, the air will be sucked back into the major 
pipe, so no dangerous release of VOC would occur. 

The air might contain a large amount of methanol (up to 2000 ppm) along with other more 
dangerous VOCs. Fortunately, methanol is highly soluble in water and could be removed by 
scrubbing. The water is sprayed on the scrubber packed bed, which is about 6 feet high and is 
filled with porous media enhancing mass transfer between the liquid and gas phases. 

Downstream of the scrubber, the water will be virtually free from methanol. The air will 
contain 100% of relative humidity. As will be seen in the research part of this thesis, the 
operation of corona in humidified air is more efficient, than with dry air. 

Pulsed corona module 

The design of the corona module is very similar to scrubber. Instead of a packed bed, 
there are 12 standard corona "blocks" (4 stacks x 3 blocks). The design of a corona block is 
shown in Figure 2. 
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Air and VOC 




Figure 1. A 3D-view of the trailer-mounted pulsed corona pilot plant. 



The block consists of intermitting layers of negative and ground electrodes. Such 
configuration creates a quasi-uniform volumetric discharge with minimum dead volume. The 
sprayed water plays a multiple role. It could dissolve some additional amounts of non-soluble 
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compounds such as dimethyl sulfide or terpenes (primarily alpha-pinene) in accordance with 
Henry's law. 



Figure 2. Design of a corona block. 

But water also dissolves dangerous byproducts of VOC destruction. This fact is very 
useful for applications. Instead of converting VOCs completely into final products such as CO2 
and H2O, one can spend much less energy on creation of solvable byproducts with OH groups. 
Solubility of such molecules is higher due to polarity of OH groups. Therefore, higher removal 
rates could be achieved in scrubbing of non-soluble in water compounds. Once in the solution, 
the organic material could be treated by more efficient methods such as biofiltration. 

Droplets of water capture significant amount of ozone. Whereas ozone does not react with 
molecules of VOC in the gas phase, in the water it inunediately creates additional OH radicals 
that hunt for primary VOC, dissolved in water according to Henry's law, or for byproducts. This 
also increases the overall efGciency. The effect of enhancement should be pronounced in strong 
alkaline solutions (pH > 9), because the solubility of ozone exponentially increases with pH. 



Ground 



Electrodes- 
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Mist Separator 

It is desirable to collect and recycle all the water sprayed over the corona. The droplets 
with sizes smaller than O.l mm will be drugged away by air. Therefore, the sprays in scrubber 
and corona modules should produce droplets larger than 0.1 mm. However, too large droplets 
should be also avoided because they reduce the surface for mass transfer. In the current pilot 
plant the sprays produce droplets with diameters in the range between O.I and 1 nun. 

A stack of trays with activated carbon in the mist separator will effectively retain the 
droplets. After the corona module, the stream will contain very small amoimt of VOCs. Those 
quantities are very diflicult and expensive to treat. But the carbon will retain the rest of VOC and 
byproducts on its surface and accumulate them for some additional time. It is known that the 
surface also promotes reactions of ozone with organics. The reaction rate of ozone on the surface 
of activated carbon will be high, and less than 1 ppm of VOC will be expected to leave the pilot 
plant after the mist separator. 

Control Equipment 

The advantage of pulsed corona is in the possibility to adjust its power by changing the 
pulse repetition rate (PRR). As will be described later in this thesis, the production rate of firee 
radicals in corona is proportional to the corona power. It was found by the modeling of corona 
chemistry that VOC react with free radicals according to the first-order reaction. Therefore, there 
is an almost linear dependence of VOC destruction rate on corona power. 

The control network of the corona pilot plant consists of several instruments. The Total 
Hydrocarbon Analyzer (THC, TEI Inc) is a robust industrial instrument capable of instant 
monitoring of the hydrocarbon concentration in the gas stream. The sample of air is burnt in the 
hydrogen flame in the Flame Ionization Detector (FID) that produces the electrical signal 
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proportional to hydrocarbon concentration. The signal is logged into the Data Acquisition System 
(DAQ) based on a PC-type computer. DAQ adjusts the output power of the corona control 
module via the pulse repetition rate. The Gas Chromatograph / Mass Spectrometer (GCMS, 
Varian Inc.) takes periodic sampling of the exhaust gas. GCMS work cycle is from I to IS 
minutes. Therefore, instant control over the exhaust is impossible. Fortunately, the VOC 
concentration in the major vent pipe changes very slow with cycle of about 30-60 minutes. So, 
the presented instrumentation will be sufficient for the pilot operation to demonstrate the 
technology and maintain safe operation. 

i.lO. Scientific Problems and Research Objectives in Development of the Wet Pulsed 
Corona Technology 

In previous paragraphs we discussed the hardware for demonstration of the pilot pulsed 
corona technology. The logic and optimization of its testing and foregoing operation demanded 
specific data on physics and chemistry of the discharge at conditions of the ''worst case scenario". 
The summarized peculiarities of this scenario are: 

1 . Concentrations of major components are out of the guaranteed operational range of corona. 

2. Poisoning of discharge with methanol. 

3. Large amount of water in the stream. 

4. Possible variation of stream temperature (operation of corona with or without sprays). 

5. Existence of the liquid-gas interface in the discharge. 

Some items in the list above would constitute a real technical problem. The information 
on others would bring new knowledge on plasma chemistry and physics of non-equilibrium 
atmospheric-pressure discharges. Let us state the relevant objectives of the research. The 
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testing and further application of the pulsed corona technology to treat VOC emissions in the 
forest products industry. 

Here we underline the tasks to fulfill this objective. 

1. Investigate the possibility of operating corona reactors at high (>1000 ppm) concentrations of 
VOC, particularly with large content of methanol, DMS and acetone. 

2. Study the effects of absolute humidity and temperature on the removal efficiency of VOC. 

3. Develop a model of corona chemistry using a sample substance (methanol). 

4. Using the model categorize the nature of observed effects (physical, chemical). 

5. Predict removal of the other VOC (provided the detailed chemical mechanism for these 
compounds is known). 

6. Develop recommendations for the operational regimes during testing of the pulsed corona 
pilot plant. 

7. Develop relevant figures of merit and compare pulsed corona technology with gliding arc. 

Following these objectives, a theoretical and experimental work on corona discharges and 
gliding arc will be presented in die next chapters. 



2. NON-EQUILIBRIUM LOW-CURRENT GLIDING ARC DISCHARGE 



2.1. Introduction 

The gliding arc is a high-pressure gas discharge with high electron temperature and density 
and has been proposed for chemical gas processing [2]. The arc (Figure 3) starts in a narrow gap 
between two or more diverging electrodes in a gas flow when the electric field in this gap reaches 
approximately 3 kV/mm in air [4]. Then the arc current increases very rapidly, and the voltage on 
the arc drops. If the gas flow is strong enough, it forces the arc to move along the diverging 
electrodes and to elongate. The growing arc demands more power to sustain itself. At the moment 
when its resistance becomes equal to the total external resistance, the discharge consumes one- 
half of the power delivered by the source. This is the maximum power that can be transferred to 
the arc from the constant-voltage power supply. Next, die length continues increasing, but the 
supplied power will be insufficient to balance the energy lost in heat transfer to the surrounding 
gas. The arc cools down and finally extinguishes. The next cycle starts immediately after the 
voltage reaches the breakdown value, usually just after the fading of the previous arc. (If the 
voltage is high enough, and the gap is very narrow, a new arc starts even before extinguishing of 
the old one [8]). A typical repetition rate of the arc is in the range from 10 Hz to 100 Hz and 
changes with the gas flow rate: the higher is the flow rate, the higher is the frequency. 

Unlike regular high-current thermal arcs where both gas and electron temperatures (Fand 
Te) could be as high as lOOOOK, the low-current (from 0.1 to 1 A) gliding arcs are believed to 
operate in non-equilibrium regimes (T < Te), and have low gas temperatures. 
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The arc voltage could be as high as several kilovolts. Typical vibrational and translational 
temperatures of the gas in a low-temperature gliding arc were measured to be respectively about 
2000 - 3000K and 800 - 2100 K [3]. Similar features are also intrinsic for microwave discharges. 
It shows that low-current gliding arcs could be definitely useful for plasma chemistry. 




Figure 3. Electrical Schematic of the DC Gliding Arc. 



25 

The physics of the gliding arc evolution is not well understood, especially the 
development of non-equilibrium conditions and stability at diminishing currents. Older models 
[9] did not account for variation of plasma conductivity and temperature during the gliding arc 
progress.This approach did not explain some recently discovered non-equilibriimi effects in the 
low-current gliding arcs. An interesting and important experimentally observed [10] feature of 
such arcs is an effect here after called "overshooting". It refers to unexpected continuation of the 
gliding arc cycle after the maximum power was achieved (Figure 4). 

Mario Sobacchi [13] investigated the possibility of removal of specific chemical 
compounds with a low-current gliding arc discharge. A typical value of energy consumption for 
this discharge was 1 kW-hr/m^, when the targeted value was 0.02 kW-hr/m^. It was noted that 
non-uniformity of the discharge imposes certain difficulties to treat the gas efficiently. 
Experiments were performed using three different external resistance values, i.e. 25, 50, and 100 
kW. The output voltage of the power supply was 10 kV, when the air flow rate ranged from 20 to 
80 SLM. 

Methanol (100 to 1000 ppm) was used as a substance to demonstrate viability of the 
glidmg arc technology. In some runs, water was sprayed along the flow direction. It was found 
that at 20 SLM of air and independently of VOC concentration, all methanol was removed from 
the stream at power consumption of 0.2 kW hrW. No organic byproducts were detected by Gas 
Chromatography in the exhaust stream. The results were interpreted as encouraging, showing the 
possibility for further reduction of power consumption. Unfortunately, NOx was detected in the 
exhaust at level of 3000 ppm, much exceeding all reasonable limits. It could be expected, 
however, that for more non-equilibrium arcs with lower currents or higher flow rates of air, the 
temperature of the arc channel will be lower, and reasonable levels of NOx will be achieved. 
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Figure 4. (a) Current - Total Power Experimental Dependence from [9]. The arc evolutes from 
the state of maximum current and minimum power to the critical point (ler, Pmaxi negative slope) 
where it should extinguish in accordance to previous models [6]. But these experimental results 
demonstrate existence of the arc even after the critical current was reached (positive slope); (b) 
Discharge Length - Voltage of the Gliding Arc in (a). 
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In this thesis a new model explaining the evolution of the low-current gliding arcs is 
presented. Below the factors giving rise to '^overshooting" will be discussed in terms of stability 
of major parameters of the gliding arc. Analytical and numerical solutions explain the general 
behavior of the low-ciurent gliding arc and are in a good agreement with experiment. 

2.2. PhenomenologY 

There are few parameters to consider in the analysis of the gliding arc. In this section the 
phenomenology of the discharge and a possible range of major parameters will be discussed. 

The characteristic time zi of arc formation can be found from the kinetic equation for the 
electron concentration 

—f-= k^n^n^ - nJXi (1) 
at 

where Ar, is the ionization constant, and and are electron and gas concentrations. The 
estimates for clean air and electric current of the arc 7 = 1 A give the value Ti » 1 ^s [S]. Within 
this time, low-resistance plasma is formed and the voltage between the electrodes falls. Usually 
the time of the arc formation is less than 0.1% of the total gliding arc cycle time. 

The arc moves with the gas flow. At small flow rates the velocities of the arc and of the 
gas are nearly equal. It was shown for relatively small power (200 W) that there was no 
difference between the absolute velocities of the gas (10 m/s) and that of the gliding arc [6]. At 
higher power levels (~ 2 kW and above) and flow rates as high as 830 to 2200 cmVs (according 
to our assessments based on data from [7], this flow rates correspond to the range of velocities 
from 1 1.5 to 31 m/s), the difference in velocities increases from 1 to 10 m/s. Authors of Ref.[8] 
showed how the relative velocity of the gliding arc changes with the discharge parameters. 
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Propagation of the gliding arc can be studied for two types of reactors: 1) Diverging-flow 
reactor (velocity of the gas decreases when the length of the arc increases, heat losses per unit 
length of the arc decrease along its path) [8]; 2) Uniform-velocity reactor (gas velocity is the 
same at every point of the reactor, heat losses per unit length are approximately constant) [9]. For 
both reactors there were developed qualitative and quantitative theories of arc propagation [8, 
12]. Electric parameters of the gliding arc were considered analytically [8, 9] and numerically [8]. 
Properties of the particular gliding arc can depend on the arrangement and presence of different 
elements in the electric circuit. Typically, gliding arc discharges have been studied with an 
external resistance R and a constant-voltage power supply (Figure 3). It is unsurprising, however, 
that the results obtained on each type of the reactor could not be easily compared with those 
obtained on a different reactor. Here after a semp with a uniform-velocity reactor (Figure 3) will 
be considered. 

A simple model of the gliding arc was developed from Ohm's law [9]. Power W per unit 
length (specific power) and energy losses per unit length (specific losses) were assumed constant. 
The analysis with these assumptions showed that the current of the arc drops to one half of the 
initial value during the cycle. Obviously, a substantial change of the current should change the 
arc temperature. So in the current analysis we shall take into consideration the dependence of the 
arc temperamre on current and the dependence of plasma properties on temperature and electric 
field. 

For any kind of gliding arc Ohm's law can be written in the form 

V^^RJ^^ (2) 

where R is the external serial resistance (Figure 3). 



Solve (2) for current: 

2R 

The '*+" solution corresponds to a quasi-steady state of the gliding arc column. In case of 
W = const the ' implies a negative differential resistance of the circuit 

"'^'^"-V^" (4) 

and leads to an unstable regime. 

As one can see from Eq.(3), the total power of the arc /» = >F / cannot exceed maximum 

power 

P 

^"""4/? (5) 
because if P > P„ax, the expression under the square root becomes negative and no real solution 
exists. 

At constant specific power, F changes only when the length of the arc column / increases. 
A critical (maximum) value of length 

° 4WR (6) 
cannot be exceeded. Any additional increment of the length would cause the arc instability. That 
is why in the case of a thermal arc, a fast transition to non-equilibrium (FENETRE) state may 
take place [S]. The arc would instantaneously and sharply extend. As a Onal result of transition, 
the gas in plasma would rapidly cool to temperatures about 1000 - 2000K.. Simultaneously, 
electric field £ should increase together with the electron temperature r«, and the discharge might 
be sustained for an additional period of time. This post-transition period could constitute a small 
portion of the total gliding arc cycle time [10]. But in applications, most of the useful plasma 
chemistry would take place only in this minor interval because only here would The low enough. 
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This effect called "explosion of length" was studied in [9]. According to the theory presented 
there, during such a transition the length / could be as large as 3 ler. Authors of [3] observed this 
phenomenon. Other researches [10] could not confinn the "explosion" effect and made a 
conclusion that the transition between equilibrium and non-equilibrium regimes is smooth. The 
lack of reproducibility could be probably explained by the influence of some additional 
parameters of the experimental setup that were not taken into account. 

Recent experiments [10], where the total power and plasma voltages were measured as 
functions of electric current and plasma length (Figure 4), demonstrated some stability of the arc 
in the "overshooting" regime. When the power of the arc reached maximum possible value, the 
arc did not extinguish as fast as it was expected. Instead, the power began smoothly decrease and 
the plasma channel was still visible for considerable time. To explain this phenomenon we take 
into account that the arc current drops during the cycle, and a substantial change in the current 
should alter the arc temperature and specific power. A new model of the gliding arc incorporating 
these features will be described in Section 2.3. 

It will be seen later that the theory with fV= const underestimated the range of stability of 
low-current arcs. But it can be used for assessments in cases with a limited set of given 
parameters such as R, Vo and /. 

The stability analysis of energy balance equation and Ohm's law is carried out in Section 
2.4. The energy balance equation and Ohm's law will be solved analytically and numerically in 
Section 2.S. In section 2.6 we shall discuss the results obtained using our model. 

23. Transient regime of the gUdine arc discharge 

It is possible to distinguish the transient regime in the low-current gliding arc evolution. 
The gas temperature T decreases and departs from electron temperature r«. The temperature- 
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sensitive mechanism of thennal ionization yields to the field-sensitive mechanism of direct 
electron impact ionization. Therefore, one should consider the presence of both mechanisms for 
modeling. The form of the governing balance equations would remain the same, but now the 
electrical conductivity must be dependent on the electric field as well. In the following analysis, 
we introduce a simple form of the field-dependent approximation, which still permits a linear 
stability analysis. 

Consider the arc as a symmetric cylindrical column. The electric field E is aligned with 
the longitudinal axis and is assumed to be uniform within every cross-section of the arc. All other 
plasma parameters depend only on the arc radius (Figure S). 




Figure S. Variables used to describe the arc channel. The radius of the arc is defmed by the ratio 
of conductivities equal to V (after Raizer [4]). 
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A heat transfer balance equation of a quasi-steady arc is derived in accordance with the 
Elenbaas - Heller model [4]: 



where a(T, E) and ^ r> are the electric and heat conductivity at temperature Tand electric field 

E. 

The second equation is Ohm's law in the form: 

The first term on the right-hand side of Eq.(8) is the voltage on the resistor R when the 
current J = f afT.EJE dS = a(T.E)ES flows through the network. Here S is the transversal cross- 
sectional area of the arc. The second term is the voltage on the plasma column. The electrical 
conductivity of plasma is taken as a following function of temperature and electric field [4]: 



where Ob is a constant for the given gas, k is Boltzmann constant, / is the ionization 
potential of the gas molecules. 

The major difficulty here is to determine the dependence of Te on T and £. For thermal 
plasmas Tc = f and the electric conductivity depends solely on the gas temperature. For cold 
plasmas Te depends on E only. The intermediate case under consideration, which is the most 
promising one for many applications, has not been studied in detail. For our semi-quantitative 
model it is reasonable to extend the equilibrium approach to more non-equilibrium cases by 
taking electron temperature Tc in a simple form [10]: 






(10) 
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where £, is a reference value of the electric field characterizing transition from thermal to direct 
electron impact ionization [1 1]. This value is about 1 kV/cm for air at standard conditions [4]. In 
the present work the case for which E < Et is considered. The formula (10) shows that if the arc is 
ahnost thermal, its electric conductivity mainly depends on the temperature. The effect of electric 
field is more significant or even prevailing in transitional and non-equilibrium regimes. 
A temperature dependence of thermal conductivity can be expressed as in [4] 

A(r,£) = ..,xp(-^). 

where ^ is a constant, and D is the energy of dissociation. 

The logarithmic sensitivities of the electrical and thermal conductivity to gas temperature 
can be found from Eq.(9) and Eq.(l 1): 

d\n<s(T.E) ^ I 

d\nT ~ IkT ' (12) 



^ 51nr ~ 2kT ' (13) 

as well as the logarithmic sensitivity of electric conductivity to electric field: 

d\na(T.E) ^ IE' 

ain£ kTE- ' (14) 

We shall use Eqs. (12), (13) and (14) for stability analysis in the next section. 

The radius of discharge To is defined as the distance from the center of the arc to the point 

at which the electric conductivity drops "e" times (Figure 3 and Ref. [4]): 

«^(ro) (15) 
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The area S = xii 'is expected to increase with the electric field when the discharge 
becomes more non-equilibrium. Let us introduce the temperature drop across the arc Ar= T(ro) - 
T(0). If the temperature behaves as a smooth function of r, it could be expanded in a Taylor series 
in the vicinity of r = 0: T(r) » T(0) + 1/2(3^7/81^)^^1^- After rearrangement: 

— -CiT)ri . 

T (16) 
Here C( T) is a smooth, slowly changing ftmction of T. 
At the same time, from Eq.(9) and Eq.(l5): 



^T 2kT(, E'] 



Combining Eqs.( 16) and (17): 



(18) 



where So is a field-independent coefficient. So, as we expected, the area S increases with the 
electric field. 

Integration of Eq. (7) (see [10]) with ^T/T taken as in Eq.(17) results in the following 
thermal balance equation: 



(19) 

This formula establishes the relationship among the thermal conductivity, electrical 
conductivity and the electric filed. Equations (18) and (19) of the new model will be used in the 
next sections for the stability analysis and simulation. 
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2.4. Stability analysis 

According to Ohm's law, the sum of voltage El on the arc and the voltage on the serial 
resistance R must not exceed the constant external voltage Kq in Eq.(8). The development of non- 
equilibrium conditions in the growing gliding arc eventually results in the discharge instability 
followed by the arc termination. When unstable, the arc cannot exist and propagate for a 
considerable time. The matter of interest here is to understand whether the arc is still stable in 
'^overshooting", or just decays. 

A stability analysis will be done for the thermal balance Eq.(19) taken in a time-dependent 

form: 

^ dT \67tkXT- 
pC,- = aiT,E)E- (20) 

where /> is a density of the gas and Cp is its specific heat. Substituting electric current 7 =of7'.£;£5 

into Eq.(8) 



K„=o-(r,£)£S^ + £/ . 
Let us consider the temperature and electric field fluctuations in the forms 



(21) 



Ar(/)=Ar,exp(n/) , ^^2) 

A£(/) = A£„exp(n/), ^^3) 
where Afo and A£o are the initial temperamre and field perturbations, and Q is the instability 
decrement. For a stable regime {Q < 0) the fluctuations will exponentially decay with time. 
Conversely, positive Q will characterize instabilities, i.e. the fluctuations growing in time. 

After substituting Eq.(22) and Eq.(23) into Eq.(9) and Eq.(ll), and then linearizing the 
result, we obtained the relative fluctuations of electric conductivity aand diermal conductivity A 
in terms of sensitivities Eqs.(12) - (14): 
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= <Jt-:=— +<Je— . (24) 



AA(0 AJ 

A r • (25) 
Now, apply expressions for t ) and £iE( t ) together with Eq.(24) and Eq.(2S) on 
Eqs.(20) and (21). Linearization and elimination of exp(Qt) yields 



n^ = a;|^(<T,-2-A,)^+(2 + ^r,)^j , 



(26) 



(27) 



where co =aE' / (p Cp T). 

Linearization of Eq.(2 1 ) gives the relationship between AE and AT : 

A£ ^ AT* 
E ~ (Tf + ! + ///„ T ' 
where /o = a 5 is the same as U in Eq.(6), and it is the length of the arc column when plasma 

voltage becomes equal to one-half of the external voltage Vo (the point of maximal arc power). 

From Eqs.(26) and (27): 

"H^'^f^r^*'*'^)- (28) 

This important result shows that the arc can be stable beyond the point of maximum 
power. Indeed, Q is always negative dl I < U and becomes positive, i.e. the system looses 
stability, at some critical length /cr > U- 



<rr+(2-t-/lr)(<Tg-Hl) 



(T,-(2 + >lr) (29) 
Previously, a stability analysis was carried out in [3], but the decrease of energy losses 
during the gliding arc evolution was not taken into consideration. There the instability developed 
inunediately after / = U, and la was equal to /o. In the experiments [5], the instability normally 
occurred at / significantly larger than U. Apparently, the new model is closer to reality for low 
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currents. Formula (29) does predict the "overshooting" effect as seen in Figure 4a: the total 

power of the arc reaches the maximum and then declines, and the arc is still stable. 

If the system is stable at / = /o, it would be natural to suppose the solution of Eqs. (19) and 

(21) exists in the vicinity this point. The expression under the square root in Eq. (3) stays positive 

because the temperature decreases, and so does W. Then the real solution should exist for / > k. 

At the same time, as we have mentioned above, the length cannot be indefinitely large because 

the voltage El of plasma is restricted by the constant external voltage Vo, and E grows together 

with /. From Eqs. (19) and (21): 

dE_^ E{cT^-2-X^) 

dl ~o-5/?[(ct£+1 + ///J(2 + >1^)+<t^(1-///j] * (30) 

Rewrite decrement Q using Eqs. (28) and (30) in the form 

toE{tTr-'^-K) 
aSR {dE/dl)i<TElo + + 0 ■ ^ ) 

Because the numerical value of (at- 2 - A.?-; is positive for our set of parameters, the arc 
would be stable and dE/dl would be positive, as expected. When / approaches /cr, dE/dl blows 
up. So, the instability develops at the moment when E is close to its maximum possible value and 
the length could not be any larger. 

This simple linear stability analysis revealed and underlined two important results: 1) The 
gliding arc remains stable during "overshooting" and 2) The instability appears only in the 
moment when the voltage drop on the serial resistance can no longer compensate the rise of the 
voltage on the gliding arc. These are the key points of understanding the experimental results [10] 
(Figure 4). 
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2.S. Numerical solution 

To justify the applicability of our approach, in this section we examine how closely and in 
what range of parameters our model matches with the experiment. 

Generally, the systems of non-linear equations such as Eqs. (18), (19) and (21) cannot be 
solved directly. A numerical solution will be discussed later. An approximate dependence can be 
obtained by taking a natural logarithm from both sides of Eq. (19). After rearrangement 
I + D(\ + E-/E;) f 16;rAr/lorM 

2it7'(l + £-/£,-) ~ [ aJSX' J ' (32) 
Numerically D« I and 1 < (l+£~ /Ef) < 2. Here the second term in the nominator on 
the left-hand side can be neglected. Due to the nature of the logarithm, the function on the right- 
hand side changes much slower than T. It can be treated as a constant value because in the model 
the temperature and electric field is not expected to change very much. 
Finally, 

"l + E-ZEf ' (33) 
where ^4 is a constant. This formula obtained from the balance Eq. (19) is very general, and 
extremely simple. It seems very useful and broadly applicable. The formula agrees with the 
experimental fact that when the current and temperature decrease, the electric field increases [10]. 

In the numerical solution, the ranges of £, T and 5 were found solving Eqs. (18), (19) and 
(21) simultaneously. The parameters were set in the following manner. The length / was varied 
from some minimum values (electrode gap) to the critical values la, where solution did not exist. 
The cross-sectional area S was chosen as in experiments as well as external resistance R. 
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Constants Go, Xo. / and D were found by fitting the experimental data [4]. The results of 
numerical solution of Eqs. (18), (19) and (21) are presented in Figure 6 - Figure 9. In Figure 6, 
the arc starts from large currents and small power (negative slope). Then it reaches maximum 
power P„ax at currents Jo. According to the earlier model of the gliding arc [6], it should become 
unstable and instantly extinguish. Conversely, our calculations with variable speciflc power of 
the arc give solutions even after this point for J < Jo up to P* (positive slope). This is in a good 
agreement with experimental results [10] (Figure 4a). 

The explanation of this "overshooting" effect could be illustrated using Figure 8 and 
Figure 9, showing calculated electric field strength and gas temperature for different lengths of 
the arc. Ohm's law restricts the total power of the arc by P^a- If the temperature remains 
constant, there will be no more power to support the arc after U (/o = Icr at constant power losses). 
Actually, when the arc approaches the critical point, the temperature decreases (Figure 9) as does 
the power losses per unit length {W), helping to maintain P„ax = IW. 

Through the relationship Eq.(lO) the model incorporates the change in ionization 
mechanism in the vicinity of the critical point. According to Eq. (33), the decreasing temperanire 
gives rise to stronger electric fields (Figure 8). When the field approaches £, this model 
incorporating a very simple dependence of Te on T Eq.(10) is no longer applicable. For the real 
gliding arc (Figure 4a) the electric field may be even stronger, so the arc is sustained farther. 
Also, the model is not suitable for strong non-equilibrium and for very low currents (case /i = 100 
kOhm, J < 0.1 A, in the experiment E > 1 kV/cm, it was assumed that £, ~ 1 kV/cm), where 
"overshooting" is very pronounced in experiment, but only partially seen in calculations. 
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Figure 6. Calculated power-current dependence of the gliding arc. The results are in a very good 
agreement with experimental results [10] (Figure 4). A stable solution exists beyond the point of 
maximum power Pmax. No physical solution exists for P < P*. The model is not applicable for 
small currents (J < 0.1 A, R = 100 kOhm). 
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Figure 7. Calculated arc voltage for different lengths of the arc. Circles denote the points of 
maximum power Pma and / = U. 
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Figure 8. Calculated electric Held for different lengths of the arc. Circles denote the points of 
maximum power Pmax and / = /g. The slope of the curves (dE/dl) increases with length. 
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Figure 9. Calculated temperature for different lengths of the arc. Circles denote the points of 
maximum power Pmax and / = U. The arc becomes more non-equilibrium as / increases. 
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Our analysis started from a few basic assumptions about the nature of the discharge and 
its properties. Also, only a few bulk parameters could be known or found from the experiment. 
Nevertheless, being a simple pbenomenological theory, the model and its numerical solution 
exhibit a behavior similar to the experiments (Figure 4). 

The linear analysis showed the possibility of the low-current gliding arc to remain stable 
even after the maximum total power P^tx was reached ("overshooting" effect). This result is in a 
good agreement with the experiment. The "overshooting" phenomenon was explained through 
the incorporation of temperature and Held sensitive properties (electric and thermal 
conductivities) into the model. The arc would gradually degrade toward the non-equilibrium. 
Driven by the falling current, the gas temperature of the arc, the electric conductivity and specific 
power iV also decrease. As a result, the total power could smoothly fall after reaching P„ax and 
the expression under the square root in the Eq. (3) (Ohm's law) would stay positive. 

Since Ohm's law limits the sum of voltages on the arc and on the serial resistance, the 
development of non-equilibrium conditions accompanied by the growth of electric field 
eventually results in discharge instability, and the arc extinguishes. The stability analysis showed 
that the instability appears at the moment when the voltage drop on the serial resistance can no 
longer compensate the growth of the voltage on the gliding arc. 

Comparison between the experimental results and numerical solutions shows that our 
model is applicable for the transition regime of gliding arc when the effects of both thermal 
ionization and direct electron impact ionization must be considered. 



3. NON-EQUILIBRIUM PULSED CORONA DISCHARGE 



3.1. Introduction 

Corona is a self-sustained electrical gas discharge that occurs only when the electric field 
is sharply non-uniform. The field near one or both the electrodes must be stronger than in the rest 
of the gap. This situation typically arises when the characteristic size r of the electrode to which 
the high voltage is applied, is much smaller than the inter-electrode distance d. If we consider a 
wire-into-cylinder configuration, the electric field in the space between the coaxial cylinders of 
radii r (internal cylinder) and R is given as a function of the radial coordinate x as: 



where K is the voltage between the cylinders. 

There is always some ionization in atmosphere due to high-energy particles coming from 
space. Existing electrons are accelerated by the electric field. They can ionize more molecules of 
the gas. New electrons are accelerated in the field, and so forth. Corona discharge occurs only 
when the value of the maximum electric field exceeds the breakdown field. For coaxial 
cylinders, the critical electric field of corona ignition in air is expressed by an empirical formula 
suggested by Peek (1929): 



(34) 




(35) 



43 



44 

where 5 is the ratio of air density to the density in standard pressiire and temperature. It is clear 
that a constraint is imposed on the maximum field, which occurs on the surface of the sharp 
electrode, that is the high field or "corona carrying" electrode. 

The formulae given above describe fairly well the physics of the corona with DC or 
slowly changing fields, such as in Gas Phase Corona Reactors (GPCR). For these applications, to 
be sustained at non-equilibrium state, the discharge can acquire only small power. Higher power 
would cause sparking and eventually lead to thermal arc formation. 

Natural power limitations of DC corona technology gave rise to development of the 
pulsed corona technology. In a pulsed corona discbarge, voltage is applied to the sharp electrode 
as a series of fast rising pulses. Incoming pulses generate a number of microdischarges 
(streamers); the faster the pulse, the more streamers are produced per unit length, increasing the 
power input per unit gas volume. Some parameters of streamers are given in Table VI. 

The width of the applied pulse must be minimized in order to avoid formation of spark 
discharge or thermalization of the streamers. For most geometries, then, the pulse rise time must 
be on the order of few nanoseconds and the duration of the pulse on the order of 100 ns. Pulsed 
corona is more suitable for the goals of the Forest Project, because all range of concentrations of 
VOC can be treated with the same discharge by simply changing the pulse repetition rate. 

To obtain experimental data, used in this work, a negative pulsed corona has been applied. 
The electron energy, in fact, depends on both the intensity of the electric field (i.e. on the voltage) 
and the mean fi%e path. For a fast rising pulse, the mean electric field is very high (since the peak 
voltage is high), allowing the electrons to gain enough energy for the discharge to take place; 
however, since the duration of each pulse is much smaller than the interval between pulses, the 
mean required power is low (e.g., it is in the order of few Watts for the setups used for these 
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experiments) even though the actual power of each single pulse is very high (about 1 MW). The 
negative polarity is predetermined by the existing equipment. 



Table VI. Typical Parameters of Microdischarges (Streamers). 



Lifetime 


1-20 ns 


Current Density 


O.l - 1 kA/cra- 


Transported Charge 


100- 1000 pCoulomb 


Diameter 


100-200 /urn 


Electron Energy 


1-20 eV 


Electron Density 


10'''-10'^cm^ 


Electric Field 


E/n=l-2(E/n)pa«hcn 


Gas Temperature 


300 K 



Usual DC corona is a uniform discharge. In contrary, the pulsed corona is highly non- 
uniform. All the discharge power is localized in streamers. To describe the physics of the 
discharge, one must consider a phenomenon of a streamer as a whole. Electromagnetic model of 
physics in streamers must be coupled with chemistry of excited species, and with internal 
particles (electrons, photons) transport. This coupling is very difficult because of enormous 
amount of computational work. Simplified models were developed, but they are still not adequate 
for the purposes of trustworthy chemistry prediction. Many parameters are known with great 
degree of uncertainty, such as electron cross-sections. In this thesis, a simpler way was found to 
represent the discharge physics and to decouple physical and chemical phenomena. 
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3.2. Review of Literature 

The topic of air pollution control by electrical discharges was subjected to intensive 
research for decades. In this Held, the plasma-chemical processes in the most challenging area of 
study. Developments in research technique (both software and hardware) for the last IS years 
significantly enhanced our knowledge on elementary processes in plasma. The research was 
stipulated by constantly increasing interest to plasma processing in industry and by new stringent 
environmental regulations demanding new pollution control techniques. 

The pulsed corona discharge is a relatively new technique. In pollution control 
applications, it competes with electron beam plasma technique. Current applications of both 
processes include the removal of major polluting constituents of NO, NO2 and SO2 encountered in 
the flue gas streams and vehicle exhausts. The acute problem of removal of VOC from the gas 
streams has emerged recently due to new environmental regulations. 

As mentioned, a science of discharges has a very long history. Electron beam (e-beam) 
technique has been investigated widely [14]. Due to its uniformity, the plasma generated by e- 
beams helps smdy processes in radiation chemistry (Beta-particles are electrons). The methods of 
modeling (primary Monte-Carlo) of e-beams are also very well developed. The flux of electrons 
yields the highest ionization rate per 100 eV from any of the plasma discharges. The 
disadvantage of such systems is their high capital cost and complexity. 

Li its turn, pulsed corona equipment is cheaper, but the discharge is not uniform. As 
mentioned previously, it consists of small microdischarges (streamers). The modeling is very 
complex and requires substantial computer resources. That is why the computational research on 
streamers has become more common only within the last ten years. 
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The DC corona is widely described in the literature. There is much less literature on 
pulsed corona. It is mostly on SOy and NOx removal from streams in power plants. Very little 
information is available on VOC chemistry in impulse discharges such as pulsed corona or DBD. 

A very extensive recent review [14] of air pollution control techniques by electrical 
discharges summarized the experimental facts about some trends in behavior of pulsed corona 
and dielectric barrier discharges. The paper demonstrates that this area is mostly descriptive. 
There is no agreement or no information on such major issues as 1) best polarity of corona 
(positive or negative); 2) dependence of discharge performance on temperature and humidity; 3) 
figures of merit (researchers use several ones, but it is difficult to bring the results of different 
authors to the common ground); and, finally, 4) the upper limits of concentration for which the 
pulsed corona technology is still applicable. This reflects the experimental fact that at higher 
concentration of pollutants the efficiency of discharge increases, but more power is required. 
Studying the recent literature, a developer or user of environmental plasma equipment would find 
very little guidance through this multitude of experimental facts. 

It seems that the plasma-chemical mechanism of iVO^ and SOi formation or removal is the 
best development is this area. Various independent experiments refined this mechanism, but 
some questions still remain. Evidently, this is a good starting point for the further research and 
development of corona technology. That is why the research presented in this thesis is 
concentrated on the chemistry of the pulsed corona discharge. 

To study the chemistry of DMS, methanol, acetone and terpenes in a pulsed corona 
discharge, one would use the A^Qr mechanism as the foundation or a staring point, and add some 
reactions relevant to particular compounds. However, there are some hurdles to overcome: 
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1. The mechanism for every VOC starts from the point when the concentration of radicals 
formed in the discharge is known. Obtaining this initial concentration requires solution of the 
streamer problem or finding some approximate theoretical of experimental approaches. A 
great attention is paid in this thesis to simplify this problem. 

2. Majority of chemical mechanism data is available for dry air. When reactions with water and 
OH radicals are added, as in the case of forest industry, some additional verification of 
mechanism might be required. Also, it is expedient to determine the sensitivity of major 
reaction channels to minimize computational work. 

3. The mechanisms for VOC are taken from atmospheric chemistry. Their applicability for 
corona was not clear. 

Therefore, the primary goal for the further literature research would be obtaining 
additional information on primary mechanisms (for example, from radiation chemistry), on 
radicals chemistry, and on chemistry of complex compounds, its theoretical and experimental 

evaluation. 

The properties of the pulsed corona discharge were studied experimentally by numerous 
authors, for example [49-SS]. In [SO] the pulsed streamer corona discharge in water was studied 
with the aid of optical methods. It was found that production of OH radicals depends non-linearly 
on the electric field. 

Different aspects of chemistry of inorganic gases, such as NOx, SO2, CO and CO2 were 
smdied experimentally in the works [49, 52-55]. The detailed coverage of theoretical and 
experimental results on the pulsed corona is given in work [58]. The data used in the present 
research were taken from the works on radiation chemistry [19, 21, 22, 23]. The detailed kinetic 
model of chemistry in the corona discharge is represented in work [59]. 
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The issues of formation and decomposition of ozone are iliuminated in monograph [48]. 
The detailed kinetic description of oxidation of SO2 and formaldehyde CH2O (major byproduct), 
and also the kinetics of active particles in the discharge, are taken from [43] and [46]. 

Several aspects of the low-temperature oxidation of sulfur-containing organic matter, 
panicularly DMS, are described in [IS, 16, 25-29, 31, 32, 34-37]. The need of using so 
significant quantity of sources is related to the fact that at the present moment kinetics of the 
oxidation of the sulfur-containing organic matter is not completely understood. The matched 
mechanism of oxidation of DMS under the standard conditions is given in work [17]. Interaction 
of ozone with VOC in the liquid phase is examined in [17, 26, 40-42, 25, 47]. 

bi the following chapters, the chemistry of elementary processes will be considered. The 
approach in determination of free radical concentration as an input to VOC chemistry mechanism 
will be outlined. 

3.3. Gas Discharge Chemistry 

The chemical processes in discharges start when electrons achieve some threshold energy 
sufficient to excite or ionize molecules in the gas. In strong electric fields, between collisions 
with neutral species electrons are capable of gaining the energy that is sufficient for multiple 
ionization events. The ionized species neutralize and may form free radicals. Very active radicals 
react with VOC to produce intermediate species. The processes with electrons occur very rapidly, 
while recombination is much slower. The lifetime of radicals is the longest process if compared 
with the timeframe of ionizatioi^and ions neutralization. Therefore, using the characteristic time 
as a parameter, all three stages of discharge chemistry can be distinguished (Table VII). 
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TABLE VU. STAGES IN DISCHARGE CHEMISTRY 



Process 


Example 


1. Electron - neutral collisions in corona streamers 


Dissociation 


N2 + e->N + N + e 


Electronic molecular excitation 


N2 + e N2' + e 


Formation of molecular ions 


N2 +e -> N2'' + e + e 


Simultaneous ionization and dissociation 


N2 + e -> N + N* +2e 


2. Neutralization and quenching 


Quenching of the excited state 


N2'+N2->N2 + N2 


Formation of a molecule from individual atoms 




with participation of the third particle 
absorbing the excess momentum 


N + N + N2-*N2 + N2 


Charge transfer 


N2" + H20->N2 + H20* 


With low-energy electrons 


O2" + e -> O2 
O2 + e O2" 


Ion-Ion recombination 


H3O" + NO2" H2O + NO +0H 


3. Chemical reactions with the radicals 


Among radicals 


0H + 0H-^H20 + 0 


Among radicals and VOC 


CH3OH + OH -> CH2OH + H2O 



Here ~*" denotes the electronic excited stale. 



The non-uniform nature of the corona discharge and existence of its streamer 
microstructure results in substantial challenge in modeling of the discharge. Typically, the 
researches model streamers to find electron temperatiire and to obtain the concentration of 
chemically active species (O, OH, N, A^", etc). This approach requires exact knowledge of 
electron transport cross-sections and leads to huge amoimt of computational work. To reduce the 
computational time, a set of approximations is used. This gives quite satisfactorily qualitative 
results, but the quantitative results may significantly deviate firom the experiment. Streamers 
usually have a tree-like structure. Accounting for this phenomenon would need more complex 
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modeling and much more computational time. Such modeling is very cumbersome for current 
scientific and engineering applications. 

A different approach is used for modeling of the electron beam plasma. As reported in 
[62], a beam of electrons (typically with a very narrow electron energy distribution) has the 
highest possible yield of radicals per unit of energy consumed in the discharge. The mechanisms 
and theory of interaction of electron beams (and other types radiation) with different substances 
is of primary interest in radiation chemistry. The yield of active species or radicals is usually 
described in terms of G-factors (Gross-factors), or the number of species of a particular kind 
produced per 100 eV of energy absorbed in the media of interest. Using G-factors, one can 
estimate the concentration of radicals produced under the irradiation of the given gas with a beam 
of electrons of known energy. It will be shown below how to utilize the concept of G-factors for 
accurate prediction of removal efHciency of VOC in pulsed corona discharges as well. 

3.4. Approach in Investigation of Phvsics and Chemistry of the Pulsed Corona Discharge 

Even though a huge amount of experimental data was available on destruction of VOC in 
different discharges, it was virtually impossible to predict chemistry at conditions diHerent from 
those described in papers. Systematic studies of dependence of VOC destruction on temperature 
and humidity were not carried out for the compounds of interest. The novelty of the current 
research is in the approach of treatment of chemistry of high-energy electrons. A common 
approach with G-factors used in radiation chemistry was accommodated for the pulsed corona 
discharge. Methanol has been chosen as a sample substance for validation because its detailed 
mechanism is well known. 
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The presented study of chemistry and physics in corona discharge was divided on several 
logical steps: 

1. Theoretical investigation of the discharge using the assembled chemical mechanism for 
species of air and for a simple VOC. hivestigation of applicability of the chosen discharge 
model. 

2. Using the results of theoretical investigation, definition of the matrix of experiments. 

3. Selection of the proper quality factors or figures of merit. 

4. Modification of the theory using the experimental results. 

3.4.1. Description of Bulk Discharge Parameters 

In the following analysis the discharge was considered in the cylindrical geometry such as 
"wire-in-tube". However, the presented approach can be modified for different geometry. 

Let us examme some bulk parameters of the experimental installation, used in the 
computer model, and some relationships between them. Let Z, be the length of the corona reactor 
with the inside radius a, and the volume Vr. Then 



(36) 



For our system L = 1200 mm, a = 1 1 mm, Vr = 45.6 cm^. 



Let us Qo for the gas flow rate in SLM. With temperature correction 




(37) 



where T is the temperature of gas. To = 293 K. 



The residence time of the gas in the reactor 



'midemx 



Q ar 



(38) 
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Specific energy input will be 




(39) 



where W is the electrical power. The specific power per unit volume is 

^^-y- (40) 

It can be seen that the residence time decreases at higher temperatures. A change in the 
concentration with the temperature was also considered in the model with the aid of standard 
procedures form the CHEMKIN library. The power of discharge' varies with the temperature. The 
values of power used for numerical modeling are based on those measured in the experiment. 

It is possible to assume that the energy of the pulse does not depend on the Pulse 
Repetition Rate (PRR). The pulse energy is the function of I) electric potential between the 
electrodes; 2) temperature of the gas in the reactor; 3) concentration of the components in the gas 
mixture. 

A good assumption on the linearity of the discharge power with PRR was confumed in 
actual experimental measurements. The next formula concludes the bulk model of the discharge. 



The coupling of bulk parameters of the discharge with its electron chemistry was 
managed through using concept of G-factors taken from the radiation chemistry. 




(41) 



54 

3.4.2. Elementary Chemistry 

Development of chemical mechanism would proceed in two phases. Phase 1 will include 
preliminaty investigation of applicability of a G-factor-based concept to description of pulsed 
corona discharges at simplified conditions (dry air, no VOC). At Phase 2, VOC and humidity will 
be added to the gas stream, increasing the complexity of the model. 

Current investigation, as mentioned above, uses G-factors to couple the bulk parameters 
and chemistry of the discharge. G-factors will be defmed as a quantity of particles of the type 
"A" of those appearing in the medium during the absorption by it of 100 eV of energy of 
radiation. G-factors depend substantially on the density of the media. The concept of 
"dosimeter" as a standard for measuring the dose of the absorbed radiation is used. The following 
relationship takes place: 

GiP), =G(P')o -^-—^^ (42) 

where G(P)xf and G(P')d are the chemical radiation yields of products P and P' formed in the 
sample and dosimeter, respectively; [P] and [P^ are the chemical yields of the two products in 
units of mole per unit volume (e.g., mol * m'^, or mol * L"') when the sample and dosimeter are 
exposed in the same radiation Held; pu and po are the densities of the sample and dosimeter; and 
Dd/Dm is the ratio of absorbed dose per unit mass in the sample and in the dosimeter. The ratio 
Dd/Dm is found for different types of radiation and different types of dosimeters. 

Since the data of radiation chemistry were acquired in essence for the pure gases 
(exception is nitrogen-oxygen mixtures), by the dosimeter we will imply the appropriate pure gas, 
for example, water vapor, and we will calculate G - factor for the mixture of gases, used in the 
experiments. In our case this mixture is N2 IO2 IH2Q /VOC. 
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The ratio of the radiation dose absorbed by pure substance to the dose absorbed by the 

mixture can be evaluated as 



— = Q), ^— — 



(43) 



This is the electron fraction of component 1 in the mixture, where a}\ is the mass fraction of 
component 1 in the mixture. Through iZ/A)\ is marked the ratio of a total quantity of electrons in 
the molecule of component 1 to its molecular weight 



4? (44) 



The given formulae reflect the concept, according to which the portion of energy, 
absorbed by the given component, is proportional to its mass fraction, and the absorption of 
radiant energy occurs in essence due to interaction with the electron shells of atoms. As a whole 
the expression 

J£^i)^ (46) 

is called the electronic fraction of this component in the mixture. 

If product P can be formed under the radiation effect both from the component A and 
from the component B in the mixtiu-e, then cumulative G-factor can be evaluated as 

GiP) = GiP),e,+G{P),e,^ (47) 
where G(P)a and G(P)b are the radiation-chemical yields of P from ^4 and S and Sa and sb the 
electron fractions of ^4 and B, respectively. 
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Some authors [19] suppose the electrons with the relatively low energy (such as corona 
electrons) could effectively interact only with the electrons of the outer shells of the atoms. Only 
they should be accounted for in calculation of (Z/A) because the energy of incident particles is not 
enough to liberate electrons from the inner shells of atoms. Therefore, the contribution of internal 
electrons is practically negligible. They do not absorb the incoming radiation. 

If all the energy of the electric field would go to high-energy electrons and all these 
electrons would have the energy above a certain threshold value, one could use data for P-rays 
from radiation chemistry to calculate the quantity of radicals. However, not all the energy of the 
electric field goes to high energy electrons. The distribution of electrons is far from narrow with 
a major fraction of low energy electrons. The electron distribution function must be calculated 
using Boltzmann equation. 

In this thesis it is proposed to find the fraction of energy that goes to high energy electrons 
by measurements of concentration of long-living chemical species. This approach requires having 
a very well known detailed chemical mechanism with minimum intermediate steps to minimize 
inevitable errors. 

In our simple kinetic model the temperature is accounted for through the introduction of 
parameter a. which represents a fraction of the discharge power consumed in formation of active 
species (Figure 10). It must be the same for all the processes at given temperature. This fact was 
verified experimentally. At each pulse in corona there forms a specific quantity of primary active 
particles O, 0\ N, N', H, OH: 



57 

where gP" is the G-factor of /-th particle, calculated for this mixture compositioii, Sunp is the 
energy of pulse, equal to the ratio of power W to the Pulse Repetition Rate Q-. 

W 



X = measured discharge power Production of species with high-energy electrons 

^ . can be estimated using G-factors for e-beams 




high energy electrons. Only they participate in 
chemistry. 



Figure 10. Concept of a parameter. 



The increment in concentration of the particles of the i-th kind after pulse 

fir* ^^^^ 

When the pulse ends, a certain number of neutral particles are considered to be present in 
the discharge. The model describes the chemistry of these particles in humid air with respect to 
VOC. The algorithmic realization of model is based on integration of the Cauchy problem with 
the initial conditions determined by the species concentrations at the moment after the pulse on 
the time interval equal to the gap between the adjacent pulses. 
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In Phase I of chemistry investigation, a mechanism of NOx formation was chosen for 
determination of a(T) (see Table VIII). 



TABLE VIII. REACTION MECHANISM OF NOx FORMATION TO DETERMINE a(D 



Processes with Oxygen 


O: + e -> O2' + e (G = 2.0) | O: + e -> O2' +2 e (G = 2.07 [22]) 


O2 + e -» 0 + O' +2 e ( G = 1.23. In this case the fraction 

O2 + e -> 0* + 0* + 2e of each of the channels comprises approximately SO %.) 


02'->0 + 0' 


O' + O2 -> O2* + 0 


O2' + 2O2 -* 0/ + O2 


04' + 02"->202 + 20 


0 + 0j + M-*03 + M 


0' + 02 + M->03 + M 


Total G (O + 0' ) = 2G( O-.* )+4G( O* ) + 2G( 0-.' ) = 13.06. The fraction of excited atoms O' is about 40 

%. 


Processes with Nitrogen 


N2 + e-»N2* + e + e(G = 2.0) 


N2 + e ^ N2' + e 


N2 + e-»N + N'+e (G = 2.36) 


N2 + e->N' + N" + e 


N2 + e->N(N') + N* (G = 0.7) 


N* + 2N2->N3* + N2 


N2' + 2N2-»N4* + N2 


N3* + e-»N + N2' 


N4' + e-»2N + N2 


N + N + M->N2 + M 


Processes with Oxygen and Nitrogen 


+ O2 -* NO* + O' 


NO + NO + 0: ^ NO, + NO, 


N* + O2 O2* + N 


NO2 + O2 -+ NO + O3 


N2' + O2 -» O2* + N2 


0 + NO2 + M -> NOj + M 


N4* + 02->02* + 2N2 


O + NO2 -»N0 + 02 


0'+N2->0 + N2 


O + NOj-^Oi + NOi 


N' + 02-*N0 + 0 


NO + N03-»N02 + N02 


0 + N2 ^ NO + N 


NO2 + NO3 -» NO + NO2 + O2 


N + O2 ^ NO + 0 


NO2 + NO3 + M N2OJ + M 


N' + NO -> N2 + 0 


N205->N02 + N03 


N' + N02 -»N20 + 0 


N + Oj-»NO + Oi 


N' + NiO-^NO + N, 


NO + 03->N02 + 02 


N" + NO2 ^ NO + NO 


NO2 + 03->N03+02 


N + N0-»N2 + 0 


0' + N20->N0 + N0 


N + N02->N20 + 0 


0 + N0 + N2->N02 + N2 


N+N02->N0 + N0 


0 + NO + 02->N02 + 02 
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The yield of nitrogen oxides depends on the quantity of excited nitrogen atoms generated 
in the discharge. Thus, the a(T) dependence could be reconstructed from the experimental 
temperature-concentration data from corona operation in dry air using a successive iterations 
method. Determination of alpha for the particular experimental setup will be given in the 
"Experimental Results" section of this chapter. 

The next phase of chemistry investigation (Phase 2) assumes more realistic situation when 
VOC and humidity are both present in the stream. In addition to oxygen and nitrogen primary 
reaction channels, in the presence of water there will be additional reactions associated with OH 
radical formation from HrO molecules and interaction of OH radicals with other active species 
(see Table IX). The relevant values of G-factor in the pure water vapor are 3.2 {OH), 0.6 (//). 
0.65 (O) and 0.06 (Hj) [19, 22]. 

For purposes of the analysis of experimental data and simulation of the process of VOC 
decomposition in the pulsed corona discharge, one should know the detailed mechanisms of 
oxidation processes of the VOC. This process is initiated and takes place with the participation of 
the active particles made by the corona. The primary distinctive features of such processes are the 
conditions at which they occur, i.e. pressure and temperature. In our case the pressure is equal 
atmospheric, and temperature varies in the range from 20°C to 230°C. Under such conditions we 
will consider the processes of the low temperature oxidation, which are principally in a number 
of moments different from the well known processes of combustion, i.e. oxidation at a high 
temperature. 

One of the important differences in the cold oxidation from the combustion is that at low 
temperatures the action of active particles (radicals) leads to formation of metastable molecules 
and radicals, which practically do not exist at high temperatures. Thus, for instance, at low 
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temperature interaction of formaldehyde C^jO with the radical OH can occur according to the 
reaction: 

CH2O + OH -> H + HCOOH , 
which does not occur in the processes of combustion. 



TABLE IX. REACTIONS WITH O// RADICALS AND WATER MOLECULES 



H2O + e -> H2O* + e + e 


H2O + e -> QIC + H + e + e 


H2O + e->ir + OH + e + e 


H2O + e->H2* + 0 + e+ e 


H2O + e-^H2 + 0* + e + e 


H2O* + H2O H3O* + OH 


Oir + H20->H30" + 0 


H" + 2H20-»H30" + H20 


0* + H20->H20* + 0 


H + H + H20->H2 + H20 


H30"(H20)„ + H2O + M H30"(H20)„+, + M 


N2" + H20 ->N2 + H20^ 


0' + H2O -> OH + OH 


0' + H2O 0 + H2O 


0' + O2 0 + O2 


0' + N2 0 + N2 


0' + O2 + M O3 + M 


0' + O3 O2 + O2 


0' + O3 O2 + 0 + 0 


0H + 0H-»H20 + 0 


OH + OH + M H2O2 + M 


OH + H2O2 H2O + HO2 


OH + HO2 H2O + O2 


0 + OH -> H + O2 


OH + O3 ^ HO2 + O2 


HO2 + 0 OH + O2 


0 (O') + H2O2 HO2 -t- OH 


OH + O3 HO2 + O2 


CH3O + O2 -► CH2O + HO2 


HCO + HO2 CO2 + OH + H 


CH2OH + HO2 -> CH2O + H2O2 


HO2 + HO2 + M H2O2 + O2 + M 


HO2 + O3 OH + O2 + O2 


H02 + 0->OH + 02 


H + OH + M H2O +M 


H + 0H-^0 + H2 


OH + H2 ^ H2O + H 


H + O3 OH + O2 


H + H2O2 H2O + OH 


H + H02-)>0H + 0H 


H + HO2 ^ H2 + O2 


H + HO2 H2O + 0 


H + O2 + M -> HO2 + M 


0' + H2O H2 + O2 


N0 + 0H + M->HN02 + M 


OH + HNO2 ^ NO2 + H2O 


OH + HNO-»H20+NO 


OH + NO2 + N2 HNO3 + N2 


OH + NO2 + 02^ HNO3 + O2 


OH + NO2 + H2O -> HNO3 + H2O 


OH + HNO3 NO3 + H2O 


0H + N20->HN0 + N0 


OH + N-*NO + H 


N0 + H02->N02 + 0H 


N2O5 + H2O ^ HNO3 + HNO3 
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Furthermore, the temperature dependence of the reaction rate constants at high 
temperature regions cannot be directly extrapolated in the region of low temperatures (even for 
the non-threshold reactions with the participation of radicals). 

Chemistry of cold oxidation is studied comparatively worse than chemistry of 
combustion. However, in recent years an essential progress was made in this area, stipulated, first 
of all, by the interest of researchers to the processes of atmospheric chemistry involving different 
organic and inorganic compounds and leading to the greenhouse effect, destruction of ozone 
layer, photochemical smog, acidic rains and so forth. The data of atmospheric chemistry can be 
used for simulation of VOC decomposition in non-thermal discharges. However, one must 
consider several limitations of this approach. 

In numerical and experimental studies of chemical processes in atmosphere, the 
concentration of active particles (radicals) is usually assumed constant and equal to some average 
values. Those values depend on the rate of atmospheric processes (mainly photolysis) and on 
ambient temperature. It is frequently assumed that chemistry of VOC oxidation, to a certain 
extent, can be isolated from chemistry of active particles in the atmosphere. When modeling 
VOC oxidation in the discharge, the formation of radicals under the electron irradiation can be 
also examined separately. However, some reactions of radical losses with each other must be 
accounted for along with the reactions of radicals with VOC. That is why die reaction of OH 
radicals with each other 

OH + OH-^HjO + O 

k 2 * lO'^ mol/ cm^ sec, must be included in the mechanisms of VOC oxidation in the 
discharge. 
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There is a practice in atmospheric chemistry not to consider some active particles because 
of their short life. In contrary, in the discharge, due to large specific energy inputs the role of 
these particles is significant. Ozone formation, for example, starts from the excited state of 
atomic oxygen. The reaction mechanisms used in atmospheric chemistry account for this state 
only indirectly through the concentration of secondary active species (ozone and OH). In the 
discharge it is possible to have direct interaction of atomic oxygen with VOC molecules. 

Mechanism of oxidation of some volatile compounds in the discharges was studied by 
several authors. Kushner [43, 44] investigated the mechanisms of oxidation of formaldehyde 
CH2O and SO2 in the discharges. The rate constants of the corresponding reactions were in a 
good agreement with the data from atmospheric chemistry. At the same time, detailed 
mechanisms of oxidation of more complex substances, such as dimethyl sulfide, a-pinene and 
acetone, are presently unknown. Here the usage of atmospheric chemistry mechanisms may be 
justified. The average VOC concentration in the atmosphere is three-order smaller than in our 
system. 

The data collected by the group of atmospheric chemistry from University of Leeds (UK) 
[56, 57] and the group of Saltelli [17] were used for modeling (Figure 11). The following 
mechanism of the low temperature oxidation of methanol included the data from atmospheric 
chemistry and on the reactions of radicals HCO, HCOOH, CH2OH, CH3O with O, OH and O2 
(Table X). The full mechanism is given in the Appendix. 



TABLE X. LOW-TEMPERATURE OXIDATION MECHANISM OF METHANOL 



CH3OH + OH ^ CH2OH +H2O 


CH3OH + 0 CH2OH +OH 


CH3OH + OH ^ CH3O +H2O 


CH3OH +0 ^ CH3O +OH 


CH2OH + O2 -> CH2O + HO2 


CH3O + 02-* CH2O + HO2 


CHzO + OH-^HCO + HaO 


CH2O + OH HCOOH + H 


HCOOH + OH -> CO2 + H2O + H 


HCO + 02->CO + H02 


HCO +0 ^ CO2 + H 


HCO +0 CO + OH 


HCO +OH CO + H2O 


HCO + HCO CH2O + CO 


HCO + HCO H2 + CO + CO 


C0 + 0H-»C02+H 


C0 + 0 + M->C02 + M 
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Figure 1 L General mechanism of VOC oxidation. 
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3.43. Experimental Matrix 

The first version of assembled mechanism of methanol oxidation (in the Thesis only the 
final, tested mechanism is given) was used for preliminary calculations. They showed that the 
given set of chemical reactions not always satisfactorily reproduces earlier experimental data 
[13]. It was suspected that atmospheric reaction rates are not adequate for plasma discharges. In 
this case it would be necessary to localize at least a part of the mechanism ("reaction channel") to 
which the destruction is the most sensitive at particular conditions. The discharge physics, in its 
turn, might also be the case. Both modeling and experimental results must be available to make a 
distinction. To minimize the uncertainty, some parts of mechanisms must be "switched off' in 
artificially created conditions such as pure nitrogen, dry air, or humid nitrogen. Table XI defines 
the objectives for each set of experiments. The pulsed corona discharge would be applied to each 
of the gas compositions. 



TABLE XI. EXPERIMENTAL MATRIX 



Gas Composition in tiie Pulsed Corona Reactor 


To be Found 


Dry Air 


NOx concentration, alpha (T) 


Dry Air + VOC 


Role of "O" reaction channels 


N2 + VOC 


Role of "N and N2+" reaction channels 
radicals 


N2 + humidity 


Role of "OH" radicals 


Dry Air + VOC + humidity 


Effect of humidity and temperature 
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The other variables for experiments, in addition to gas composition, were the temperature 
and corona energy. The latter could be adjusted through changing the Pulse Repetition Rate. In 
the experiments with humidity it was varied to simulate the Forest Project conditions. Table XII 
shows relationship between the experimental conditions in air. The experiments would be carried 
out at five basic temperatures: 40, 80, 130, 180 and 230**C. For the lowest temperature (40''C) the 
absolute humidity could not exceed the saturation condition (about 2 abs %). In nitrogen a fewer 
experiments with different humidity are thought to be sufficient (0 and 2% of absolute humidity). 



TABLE XII. EXPERIMENTS: HUMIDITY - TEMPERATURE MATRIX FOR AIR 



Pulse 
Repetition 
Rate, Hz 


Temperature, "C 


40 


80 


130 


180 


230 


100 


2% abs 
humidity 


0,1,2,5,10% 
abs humidity 


0,1,2,5,10% 
abs humidity 


0,1,2,5,10% abs 
humidity 


0,1,2,5,10% abs 
humidity 


500 


2% abs 
humidity 


0, 1,2,5, 10 
% abs 
humidity 


0,1,2,5,10% 
abs humidity 


0,1,2,5,10% abs 
humidity 


0,1,2,5,10% abs 
humidity 


1000 


2% abs 
humidity 


0,1,2,5,10% 
abs humidity 


0,1,2,5,10% 
abs humidity 


0,1,2,5,10% abs 
humidity 


0,1,2,5,10% abs 
humidity 



Three basic PRR (100, 500 and 1000 Hz) were chosen to reflect two possible cases of 
discharge formation. When PRR is low (100 Hz), the new discharge does not "feel" the trace 
from the previous impulse because there are almost no active particles in the discharge gap. The 
typical time for the discharge traces to disappear is about 1 ms. Therefore, at 500 Hz and 1000 
Hz, the previous discharge would "help" the next discharge to form. Very often this causes a 
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premature breakdown. The corona becomes a spark, or may even ignite an arc (at very high 
PRR). All this must be considered and carefully monitored. The PRR of 500 Hz represents an 
intermediate case to help see a trend in transition from 100 to 1000 Hz. 

3.4.4. Figures of Merit 

There are too many parameters in the corona discharge that might influence the removal 
efficiency. Major of them are: power per impulse, voltage, pulse repetition rate, temperamre, 
humidity, residence time, concentration of WOC, composition of the carrier gas. The recent 
survey [14] on plasma echnologies available for abatement of volatile organic compounds 
showed the results of different authors in terms of several quality factors. 

The discharge parameters may be accounted for in many different ways. It is believed that 
quality factors, or figures of merit, including these parameters, must reflect the following 
experimental and theoretical results: 

1 . The Destruction and Removal Efficiency (DRE) exponentially depends on temperamre. 

2. During the detailed analysis of corona and by comparing the e^erimental data, we found that 
methanol is destroyed though die first-order chemical reaction with respect to medianol and 
radicals. 

3. Lower concentrations of VOC require higher energy inputs per unit of volume. The larger the 
concentration, the less energy is required to destroy every ppm of VOC. 

A comparison of commonly used quality factors is given in Table XIII. 
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TABLE XIII. FIGURES OF MERIT FOR VOC DESTRUCTION 



Name 


Formula 


Dimension 


Comment 


Destruction and 
Removal 
Efficiency (DRE) 


Q-c 

Q 




Does not account for 
consumed energy 


Specific Energy 
Input (SEI) 


Power* (Residence 
Time)/(Flow Rate) 


Watt*hour/m^ 


Does not account for 
concentration 


Energy Cost 


SEI 
DREC„ 


Watt*hour/m^/ppm 


Does not account for 
exponential dependence of 
DRE on power 



The parameter called "energy cost" seems to be the best to reflect the energy cost per ppm 
of VOC. But it does not show the actual percentage of the destroyed VOC. The matter of interest 
for practical applications is the degree of destruction of 90% and above. In order to find how 
much energy is needed to destroy 90% of initial quantity of VOC, one should make numerous 
experiments where he would vary the corona power. Finally, at certain power the DRE of 90% is 
achieved and energy cost is recorded. 

This long and inefficient approach was used in several works to establish DRE-power 
dependence. The amount of experimental work can be drastically reduced if the energy cost can 
be predicted from the experimental data for an arbitrary power. For the purposes of forecasting, a 
dependence of DRE on power and temperature is needed. 

The new approach is proposed in this thesis. Preliminary investigation showed that VOCs 
(in our case it was methanol) are destroyed though the first-order chemical reaction with respect 
to VOC and radicals. Then the bulk chemical kinetics equation can be written in the form: 
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^ = -k,C^JpH]-k,C^^^]-k,C^^^ = -mC^) cciT) C^W (51) 

Here C is the concentration of VOC, ki are rate constants for radicals, W is the corona 
power, P is the concentration-dependent coefficient. The solution of this equation can be 
immediately found: 

C„«(O = C_(O)exp(-)0(r,C_)a(r)fFr) , (52) 

where t is the residence time. 

Destruction and removal efficiency is then obtained: 

DRE^X- exp(- PiXX^oc) cc(T) Wt), (53) 
Using the last equation, it is possible to estimate the energy cost at DRE = 90%. The 
obtained assessment for the power will be called Projected Energy Cost (PEC), because the value 
of power for a given DRE (not equal to 90 %) is projected (predicted) for DRE = 90%. The value 
of PEC will be defined as following: 

SEI ln(l-0.9) 

PEC = 2^ (54) 

0.9 C„ In(l-D^) ^^^^ 

Projected Energy Cost has the same dimension as regular Energy Cost (W*hr/ppm/m^). 
At DRE = 90% they are equal. 
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3.S. Eiperimentai Setup 

The corona discbarge was generated in a wire- into-cy Under coaxial electrode system. The 
inner, high voltage electrode is a 0.5 nun diameter InconeP wire. The Inconet* alloy has been 
selected for its superior heat and oxidation resistance. The outer electrode is a 1.2m long 
cylindrical stainless steel tube with a 22.2 mm internal diameter. The iimer surface is polished to 
avoid localized enhancements of the electric field due to roughness of the electrode and to 
provide a uniformly distributed discharge. A tape heater created a uniform temperature Held in 
the reactor. The inner wire is held and centered by means of two Teflon® holders, located at both 
the ends of the external electrode. These pieces serve for hermetic sealing and for insulation of 
the discharge region. K-type thermocouples are attached to both inlet and outlet connections of 
the reactor, with the purpose of providing temperature data. Power is supplied to the plasma 
source by a thyratron-based power supply, with 100 ld2 internal resistance and no-load voltage 
values from 0 to 20 kV. Pulses of about 100 ns duration and 10 ns rise time are applied; they are 
transmitted to the central wire electrode by means of a SO Q high-voltage, coaxial cable. The 
power of the plasma source is varied from 0.4 to 10 W by controlling the high voltage pulse 
amplitude (typically 1 1 kV) and pulse repetition rate. Voltage is measured using a high voltage 
probe (model P6015A lOOOX, Tektronix Inc.), current measurements are performed using a 
standard current monitor (Current monitor #41 1, Pearson Inc.). 

In experiments, temperature, humidity and pulse repetition rates were set according to 
the matrix of experiments. Temperature values were 40, 80, 130, 180 and 230°C. The humidity 
was 0, 1, 2, S and 10 absolute percents. Standard pulse repetition ates were 100, 500 and 1000 
Hz. The gas flow rate through the corona reactor was 2 SLM, resulting in a residence time of 
about 13 seconds. 
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An overall schematic of the experimental setup is shown in Figure 12. Typical VOC 
compositions of tested streams were simulated in the laboratory by blending a combination of 
gaseous and liquid species to generate industrial- like streams. Inlet gas flow rate was metered by 
MKS precision mass flow controllers. The air for the experiments is provided by a dry air 
supply (dew point -70°C), that ensures the removal of residuals of water and organic compounds, 
which could be present in the air line. Liquid flows of water and VOCs were supplied to the 
system by two HPLC Pumps (models I and II, Lab Alliance Inc.) with an operating range from 
0.0 1 to 10 mi/min and 2% accuracy. The pumps use step motor driven pi.stons, ensuring good 
accuracy and repeatability, bi order to evaporate the liquid components and to properly mix the 
different chemicals, a stainless steel mixing chamber, heated up to 160°C, was placed before the 
reactors. All stainless steel connecting lines were heated by means of heating tapes to prevent 
condensation of reactants and products. The developed supply system provides reliable control 
on the chemical composition of the process streams and allows for varying of VOC concentration 
in a fairly wide range. 

Gas composition was determined by a gas chromatograph (Varian Inc, Model 3600), 
analyzing input and output samples. For the sampling, a part of the gas stream from the reactors 
was directed to the GC, collected in a 0.2S ml gas-sampling loop, injected into the capillary 
column (Supelcowax-10) by means of compressed air operated valves, and analyzed with a 
Flame Ionization Detector (FID). NOx evaluation was done by a chemiluminescence analyzer 
(Thermoenvironmental Instrument 42H). Certified gas mixtures were used for the calibration of 
the NOx aiudyzer and the gas chromatograph. Electrical characteristics were acquired with a 
digital oscilloscope (HP 546I26B, 500 MHz, 2GSa). The digitized waveforms were then 
transmitted to a computer and analyzed with the HP 34801A-BenchLink Scope software. 
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Vacuum 




Figure 12. Experimental setup. 
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3.6. Experimental Results 

3.6.1. NOx measurements 

The measurements of NOx in dry air were needed to reconstruct the dependence a(T), 
Figure 13 shows the experimental values of total NCh-NO: concentration at PRR of 100, 500 and 
1000 Hz. The function a=f(T) should be the same for all PRRs if there is a linear dependence of 
NOx concentration on corona power . As you can see from the experiment (Figure 14), this 
conditions holds for the given corona reactor. The maximum pulse voltage was 1 1 kV. 

For simulation, a(T) is taken in the form of a piecewise linear function. Its nodal points 
are shown in Table XIV. The values of are also plotted in Figure 15. 



TABLE XIV. NODAL POINTS OF a(T). 





80 


130 


150 


180 


210 


230 


a(T) 


O.l 


0.12 


0.19 


0.35 


0.61 


0.72 



In the next chapter we shall compare the experimental temperature dependence of 
methanol removal efficiency with the results predicted by our theoretical model using 
experimental a(T). Theory predicts that the increase of a with temperature is due to coupling of 
energy of electrons and electric field (E/n). Here "n" is the concentration of molecules in the gas 
phase per unit volume. 
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Figure 13. Dependence of NOx concentration on temperature at different PRR. 
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Figure 14. Dependence of NOx concentration on Pulse Repetition Rate at different temperatures. 



74 




Figure IS. Dependence of a on temperature. 

3.6.2. Experiments with methanol 

The initial concentration of methanol was 1000 ppm in majority of experiments. For dry 
air and nitrogen experiments, the concentration of methanol was 2740 ppm. It could not be lower 
due to limitations of the available experimental equipment. The matrix of experiments was 
executed, but not all experiments were possible to complete due to premature arcing at higher 
temperatures, PRR or VOC concentration. In all the experiments the maximum pulse voltage was 
11 kV. 

A pulsed corona discharge in air produces very little light, which is barely seen with a 
naked eye even in complete darkness. Luminosity of the discharge increases with PRR, but not 
with temperature. In contrary, in nitrogen (pure or with water), the discharge is very bright, and 
more discharge energy is converted to heat. Some temperature increase (up to 40°C) was 
observed at 1000 Hz. The experimental DRE-Temperature and PEC-temperatuie dependencies 
are shown in Figures from 16 to 26. 




Figure 16. Removal of methanol (2740 ppm) in dry air; a) ORE vs Temperature; b) PEC vs 
Temperature. 
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Figure 17. Removal of methanol (1000 ppm) from air with 1% of absolute humidity; a) DRE vs 
Temperature; b) PEC vs Temperature. 
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Figure 18. Removal of Methanol (1000 ppm) from air with 2% of absolute humidity; a) DRE vs 
Temperature; b) PEC vs Temperature. 
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Figure 19. Removal of Methanol (1000 ppm) from air with 5% of absolute humidity; a) DRE vs 
Temperature; b) PEC vs Temperature. 
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Figure 20. Removal of Methanol ( 1000 ppm) from air with 10% of absolute humidity; a) DRE vs 
Temperature; b) PEC vs Temperature. 




Figure 21. Removal of methanol (2740 ppm) in dry nitrogen. 




Figure 22. Removal of methanol (1000 ppm) in nitrogen with 2% abs humidity. 




Figure 24. EffTect of humidity on removal of methanol (1000 ppm) in nitrogen (1000 Hz). 
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Figure 25. Very high concentration of methanol (5480 ppm) in air a) DRE vs temperature; b) 
PEC vs temperature. 



83 




Figure 26. Effect of humidity on the removal of methanol (1000 ppm) in air at 1000 Hz. 



The difference of behavior of corona in air and in nitrogen can be explained as following. 
In air molecules of oxygen are prone to electron attachment. This leads to quenching of 
avalanches and increase of maximum allowable voltage. In nitrogen there are more excitation and 
ionization events. The active nitrogen species emit bright light. The streamers are longer, and the 
breakdown is easier to achieve. 

Some peculiarities with gas chromatograph measurements have to be highlighted. It is 
well known that the response of a FID detector is sensitive to the amount of water in the analyzed 
gas. There is also a temperature effect. Because the experiments were done at variable humidity 
and temperature, the calibration must be done at the same conditions. The calibration mixtire 
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contained no water and was at room temperature, so the calibration with the standard (mixture of 
VOC in the bottle) worked only for dry air and at room temperature. 

Fortunately, the concentration of VOC in the stream was controlled by a precision pump 
and gas flow rate was set by a mass flow controller. Therefore, the standard is given by the 
experimental setup. The measurements of concentration were done in two stages. At actual 
experimental conditions (temperature, humidity, VOC concentration) the known concentration of 
VOC was measured (10 times) at the exit of the reactor. The averaged value was taken as a 
standard. Then the measurements with corona are taken (from 10 to 15 chromatograms at each 
point), and relative change of the FID signal is calculated. These procedures gave relative error in 
determination of DRE not larger than ± 10%. 

The measurements of temperature were done by precision K-type thermocouples. The 
variation of temperature due to control adjustments did not exceed ± l^C. 

3.6.3. Experimeits with Acetone 

The experiments with acetone were needed to obtain the raw data for the following 
modeling. The gas composition corresponded to that at paper mill site, with lOORH and 
temperature about 40°C. Air contained 1000 ppm of methanol and 2% absolute humidity. This 
value of humidity was known from the previous experiment with medianol to be the optimum for 
the given corona reactor. 

The arcing happens at lower temperature and at lower PRR than in air mixture with 1000 
ppm of methanol at the same humidity. The results of experiments are shown in Figure 27. No 
traces of methanol (expected major byproduct) were found. 
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Figure 27. Experiments on removal of acetone (1000 ppm) from air streams, 2% absolute 
humidity; a) Destruction and removal efficiency versus Temperature; b) Projected Energy Cost 
versus Temperature. 
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3.6.4. Eiperiments with PMS 

The experiments were carried out in dry air. Dimethyl sulfide is barely soluble in water, 
so it was impossible to use our standard technique with evaporation of mixture. The substance 
was pumped directly into the evaporation chamber. The odor of dimethyl sulfide is unbearable, 
and even short-term exposure is dangerous for humans. All necessary safety measures were 
taken. A gas mask was worn while handling of DMS. 

The experimental curves are shown in Figure 28. Some methanol (about 100 ppm) was 
observed in the corona exhaust. Concentration of DMS, as expected, was coupled with 
concentration of methanol. At high removal rates of DMS, the amount of methanol increased, and 
vice versa. More of the discharge energy will be spend on methanol destruction. Experiments 
showed that the removal rate of DMS changes very little with temperature and even decline, 
when the energy cost increases with temperature. This is exactly opposite to what we observed in 
experiments with methanol. Therefore, to increase the efficiency of DMS removal, one should 
somehow clean the treated air from secondary methanol. A corona with a water spray will be 
feasible for such a purpose. 

As in the case with acetone, the arcing starts earlier than with methanol. Because of this 
reason, no experiments are available at temperatures above 130°C. This "cut-off' temperature is 
the same for all pulse frequencies (PRR), unlike the case with acetone. 
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Figure 28. Experiments with Dimethyl Sulfide, 1600 ppm in dry air; a) Destruction and Removal 
Efficiency versus Temperature; b) Projected Energy Cost versus Temperature. At temperatures 
above 130 °C the arcing starts. Methanol was a by-product (about 100 ppm) of DMS 
decomposition. DMS does not dissolve in water. 



3.7. Results of Modeling of Pulsed Corona Pischaige 
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The modeling was done using CHEMKIN library of FORTRAN procedures. The detailed 
reaction mechanism includes 1) NOx formation; 2) methanol oxidation; 3) DMS oxidation. Here 
NOx mechanism is the foundation for the rest of the modeling. The list of all reaction is given in 
Appendix as a printout of CHEMKIN compiler output file. 

Because the concenuration of VOC is several orders of magnitude less than the 
concentration of major constituents of humid air (nitrogen, oxygen and water), the reactions of 
VOC withelecttons are not included in the mechanism. 

From the experiments with corona in dry air without VOC, we reconstructed the values of 
a (the portion of high energy electrons). Alpha serves as a similarity parameter for all reaction 
mechanisms in the given corona setup, of course, if the experimental conditions do not change. 
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Figure 29. Relationship between NOx formation, methanol oxidation, and dimethyl sulfide 
oxidation mechanisms. 
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Figure 30. Comparison of experimental (solid lines) and simulation (dashed lines) results for 
DRE of methanol (1000 ppm) in air with 1% absolute humidity. 
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The presented mechanism of methanol oxidation includes additional reactions of radicals 
with methanol, and the following stepwise destruction of the products up to stable COz and H2O. 
At the same time, methanol is a byproduct of DMS oxidation. Therefore, to compose a DMS 
mechanism, only few reactions must be added to the methanol mechanism. Figure 29 
demonstrates the relationship between all three-reaction mechanisms. Figures 30 and 3 1 show the 
results of modeling for methanol at 1 and 2 % of absolute humidity, and Figure 32 simunarize the 
results for DMS in dry air. It can be seen that a good agreement of the theory and experiment was 
achieved. 
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3.8. Discussion 

For the following discussion, it is necessary to generalize the experimental data. The 
primary goal for the fiirther investigation is to determine the range of optimum operational 
conditions (target: minimum energy cost per ppm of VOC) for the pulsed corona pilot plant. 

In the previous chapters, the Projected Energy Cost (PEC) was proposed as a figure of 
merit to compare the experimental results obtained at different removal efficiencies. We take a 
methanol as the simplest substance to show the pattern in experimental results. 

The PEC data on methanol, obtained in several runs with different humidity, were 
arranged as a function of reactor temperature. Figure 33 demonstrates the effect of temperature. 

The same data can be plotted as a function of himiidity. The effect of humidity is 
demonstrated in Figure 34. 

Let us sununarize the major peculiarities in behavior of discharge chemistry as seen in 
Figures 33 and 34. 

1. The effect of humidity is very non-linear. At the value of humidiQr about 2 absolute percent, 
there is a minimum of energy cost for the whole set of temperatures. 

2. The temperature effect is more complex. The curves represent two different patterns: a) 
curves with smooth decrease of energy cost with temperauire (humidity of 1 and 2%); b) 
curves with a maximum of energy cost at 130 °C (humidity of 0, 5, and 10 abs %). 




Figure 33. The effect of temperature on energy cost of methanol removal (1000 ppm, 1000 Hz 
PRR) at difTerent humidities. The maximum of energy cost is observed at 130°C. There is a 
similarity of curves for no humidity and targe humidity cases (0%, 5% and 10%). 




Figure 34. Effect of humidity on removal of methanol (1000 ppm) in air at 1000 Hz. 
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In the literature we could not find any systematic data on the effect of humidity on 
removal efficiency of VOC. Majority of experiments were done in dry air (see, for example the 
work of Penetrante [62]). Methanol kinetics predicts almost linear increase of discharge 
efficiency with humidity. There should be no minima or maxima. 

The temperature effect in corona and barrier discharges were studied in several works. In 
[14], the magnitude of NO concentration in the coaxial dielectric barrier discharge decreased with 
temperature, and had minimum at about l3(fC. This result is in a very good agreement with our 
experiments. It also shows that the effect of "130°C" is not due to methanol stage of the 
mechanism. Our own analysis of humidity and temperature effects through the prism of 
chemistry in the whole range of these parameters supported this hypothesis. The temperature 
effect also seems unexplainable by means of kinetics of chemical reactions. It should be 
concluded that the abovementioned temperature and humidity effects are due to some physical 
phenomena. Here we attempt to develop physical hypotheses explaining the experimental facts. 

HvDohesis 1. Humidity effect is due to competition of two simultaneous processes . The 
minimum in the PEC-humidity curve is because of the competition of two different processes 
related to water molecules. From one side, water is needed to produce OH radicals. Chemical 
modeling shows that with the increase of water molar fraction, the destruction efficiency of VOC 
should increase. So, this is a positive role of water. From the other side, the water molecules can 
attach corona electrons. The rate of primary ionization decreases, and so does the rate of free 
radicals formation. Eventually, there will be fewer radicals to react with VOC, and the late of 
VOC removal decreases. This is the negative role of water. The balance of the abovementioned 
processes can potentially lead to maxima or minima in the ene^-humidity curve. More 
sophisticated modeling is needed to see this effect. 
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Hypothesis 2. Formation of tiny water droplets in the corona discharge. When viewing the 
temperature effect, one can see that the processes with high humidity in air (5 and 10 absolute 
percent) follow the same pattern as the process in dry air. The energy cost reaches maximum at 
130 °C in both cases. Again, this experimental fact cannot be described in the framework of the 
assembled chemical mechanism. But we attempt to explain similarities of the curves at 0, S and 
10 % of absolute humidity. 

What is the conunon ground for processes in dry air and at high humidity cases? The 
similarity would exist if there were insufficient amount of water in the gas phase. In the streams 
with high content of water, plasma can induce condensation of humidity into tiny droplets. Very 
little amount of water would then remain in the gas phase. The effect of plasma- induced 
condensation is known and studied theoretically and experimentally [63]. The presented link 
between low- humidity and high- humidity cases needs experimental confirmation with direct 
registration of not only the droplets, but also their concentration and size distribution. 

The modeling of chemistry in the discharge with the assembled methanol mechanism 
worked satisfactorily for the gas mixtures with 1 and 2 percent of absolute humidity. It is also 
feasible to extend the mechanism of methanol to DMS by adding some reactions. Therefore, 
using the method based on alpha-parameter, one should reliably predict the chemistry of different 
VOCs in the same discharge. 

The other important outcome from the corona research, in addition to understanding the 
role of temperature and humidity, is the confirmation of applicability of the pulsed corona 
technology to treatment of flue gases with high content of VOC (up to 5000 ppm). A wider 
operational range of pulsed corona reactors can be justified. For the particular case of the paper 
mill operation of the pulsed corona pilot plant, this means that the spikes of VOC concentration 
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due to the mill process cycle will be treated completely. No additional gas treatment will be 
needed. 

The buildup of methanol at high concentrations of DMS or acetone was a great concern. 
It was shown that even at 1600 ppm of DMS or at 1000 ppm of acetone, the concentration of 
methanol stays very low at about 100 ppm. This amount can be easily scrubbed by the amount of 
water that will be used for spraying in the pilot plant. If a high concentration of methanol (about 
1000 ppm) is present in the stream, it is feasible to scrub the methanol in a stand-alone scrubber, 
and then destroy insoluble compoimds in the corona reactor with a spray to wash the secondary 
methanol. 

Water can be also thought not only as a desired agent from the standpoint of process 
efficiency, but also as a necessary factor to prevent NOx formation. In dry air significant 
concentration (about 30 ppm) of NOx is formed. While destroying VOC, another undesirable 
species will be formed. The presence of some water in the gas stream helps convert these species 
to solvable HNO3. Therefore, the wet corona seems more feasible than a dry one. 

A serious difference between the gliding arc and pulsed corona is in their NOx producing 
capability. Corona produces one or two orders of magnitude less NOx than the gliding arc. 
Smaller NOx concentrations are easier to fight with addition of water. 

It has to be pointed out that from the chemical point of view our results obtained for a 
"dry" corona with 2% of absolute humidity in the stream are equivalent to the results in the "wet" 
corona with a water spray [13]. The spray does not change relative humidity whereas the air is 
already saturated with water. Because the pulse repetition rate of corona is very high and millions 
of streamers are formed per second, there is a very little probability of streamer to eiKounter a 
droplet. The chemistry with fast electrons and secondary radicals occurs in the gas phase. The 
role of droplets here is to absorb some VOC in accordance with Henry's law. 
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3.9. Conclusioas and Recommendations for Future Work 

The program of experiments aiming to enhance scientific and engineering knowledge 
required for successful testing and further application of the pulsed corona technology to treat 
VOC emissions in the forest products industry is completed. According to the tasks of the 
research objectives, the fcilowing conclusions can be made. 

Corona reactors can operate at high (>1000 ppm) concentrations of VOC, particularly 
with large content of either methanol or DMS, or acetone (each taken separately). Large amounts 
of primary methanol (about 1000 ppm) must be removed in a stand-alone scmbber in the very 
beginning of the abatement process. The amount of water that will be used for spraying in the 
corona reactor is sufficient for removing of secondary methanol that was formed during oxidation 
of DMS and acetone. 

A Projected Energy Cost (PEC) was proposed as the figure of merit to compare results on 
removal of VOC at different gas compositions and regimes of the pulsed corona discharge. 

Temperature and humidity have profound effect on DI^ of VOC. Destraction and 
Removal Efficiency increases with temperature, but there is a maximum of PEC at 130°C. 
Niinimum PEC is observed at 2% absolute humidity. It is economically feasible to humidify 
contaminated air streams for the following pulsed corona processing and keep the temperature 
either below 80°C or above 200°C. A substantial role of water is recognized to prevent NOx 
formation. 

A plasma-chemical kinetic model is developed to predict removal rates of methanol and 
DMS in pulsed corona reactors. A fraction of high energy electrons (a) was used as a similarity 
parameter for primary processes in the discharge chemistry. Together with Gfactor data from 
radiation chemistry, alpha provided necessary linking of the bulk plasma parameters and 
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concentration of free radicals. A good match of experimental results and calculation was 
achieved for methanol and DMS. 

It was found that the observed dependence of PEC on humidity and temperature could not 
be described in terms of only chemical kinetics. Two physical hypotheses were proposed to 
explain these phenomena. 

For the future work, the experiments on the industriatscale or pilot pulsed corona reactors 
are essential. An experimental program similar to that presented in the thesis should be fulfilled. 
This will generalize or limit the application of the theory of similarity parameter a in linking bulk 
plasma parameters and chemistry. The elementary chemical mechanisms for all the substances of 
interest in the Forest project should be assembled and verified with experimental data. 

The proposed hypotheses about the form of dependence of PEC on temperature and 
humidity must be supported with additional modeling of streamers and experimental results on 
quantity and size distribution of the droplets which formation is suspected in the pulsed corona 
plasma with high absolute humidity. 



4. CONCLUSIONS 



The gliding arc and pulsed corona plasma technologies for abatement of Volatile Organic 
Compounds in the forest industry were considered theoretically and experimentally. 

1 . The effect of "overshooting" in transitional regimes of low>current gliding arc was 
explained in terms of discharge stability. 

2. A concept of the mobile pulsed corona pilot plant was described. 

3. A program of experiments to investigate the operation of the plant at conditions close to 
the existing in the paper mill streams was proposed and implemented. 

4. A plasma-chemical model of methanol and dimethyl sulfide destruction in the pulsed 
corona discharge was developed. 

5. Projected Energy Cost value was proposed as a figure of merit to characterize the 
efficiency of the corona reactor. 

6. Optimum operational regimes of corona reactor operation were found with respect to 
temperature and humidity. 
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Mnta: Raactiona in Humid Air 
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